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CHAPTER I 


INTRODUCTION 


The new generation of aix-axle locomotives introduced 
In the last fifteen years has been attractive from many points 
of view, such as tractive effort, horsepower, costs, and 
ease of maintenance. However, under some operating conditions. 
In particular where the track routes are difficult or where 
the track strength conditions are low, some railroads expe- 
rience difficulties In locomotive ride, track damage, and 
occasionally, derailments. 

A serious problem developed when one particular design 
of slx"axle locomotive became widely used on passenger trains. 
These normally run faster on curves than do regular freight 
trains. After nearly five years of operations and a total 
of 21 derallments(2)*. It became evident that the problem 
deserved special attention. Various tests were conducted 
in order to find a specific reason for the derailments, but 
no single cause was identified. 


Recently, in an effort to provide a better understand- 
ing of the problem, and to develop locomotive acceptance 
procedures, extensive tests were run on specially prepared 
sections of track having geometry defects. (1) These tests 
were helpful in providing better Information on typical 
levels of wheel- rail forces and critical levels of severity 
of track geometry defects. However, such methods are costly. 

As well, the test methods provide no mechanism to predict behaviour 


* Numbers in brackets Indicate references on Pages 25 and 26. 


under operating conditions that are different from those 
of the tost. It is obvious that some testing is required, 
but the use of analytical methods to support such testing 
will greatly improve the value of the results. 

The methodology presented in this report provides 
a user-oriented analytical tool to predict operational 
safety of six-axle locomotives. 

Experimentally verified characteristics of Individual 
suspension elements, such as those obtained by Martin Marietta 
Corporation for locomotive trucks, are essential to the faithful 
representation of a locomotive. But the methodology is more 
than a mathematical locomotive representation; much informa- 
tion published by others is used for the characterization of 
several aspects such as track, wheel- rail Interaction, loco- 
motive characteristics, and derailment mechanisms. 


CHAPTER 2 


DESCRIPTION OF THE METHODOLOGY 


2.1 GENERAL OUTLINE 


Measures for assessment of operational safety are 
discussed in reference 3. Briefly, derailment can occur 
due either to flange climbing or to track gage spreading 
(vehicle roll-over does not normally occur). 

Therefore a methodology to evaluate safety of a 
specific vehicle must be concerned with wheel-rail forces. 
Since wheel-rail forces are a result of dynamic behavior, 
the following three categories of parameters must be con- 
sidered : 

I) Locomotive parameters, describing the locomotive 
mass and suspension characteristics, 

II) Operational parameters, describing miscellaneous 
conditions such as speed, tractive effort, 
lateral coupler forces, as well as track curva- 
ture and superelevation, and 

ill) Track geometry parameters describing deviations 
or defects from the nominal conditions with 
respect to cross-level, alignment, surface, 
or gage. 

These parameters are used in a mathematical model to 
simulate the response of the locomotive to the specified 
conditions. The model Is non-linear and therefore the 
solution consists of a time history for each of the system 
va r iables . 



The safety of operation is assessed by monitoring 
selected variables in the model represenMng tendency to 
derailment . 

The concept of the methodology is shown in Figure 1. 
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FIGURE 1 - METHODOLOGY OVERVIEW 







From Che current state of knowledge, flange climbing, 
gag© spreading, and permuneuC crack damage can be monitored 
through the use of the following three measures of safety s 

i) Single wheel L/V ratio representing the tendency 
for flange climbing, 

ii) Net wheelset L/V ratio representing the tendency 
for lateral track panel shift, 

ill) Truck side L/V ratio representing the tendency 
for rail rollover and gage widening. 

In order to determine the time histories of these 
measures, a software package was developed to simulate the 
dynamic response of the locomotive. Monitoring of the 
measures of safety is performed automatically by the program 
in order to determine the maximum or peak values occurring 
during a simulation. Peak values of all three safety 
measures are determined for each of the 12 wheels, each of 
the 6 wheelsets, and each of the 4 ''truck sides", A sum- 
mary of these peak values is reported after completion of 
the simulation in a tabular form. 

Options are provided for printing all important system 
variables at any instant for analysis of forces and dis- 
placements in the suspension, and for plotting time hlsr-jrles 
of any of these system variables. 
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2.2 INPUT PARAMETERS 


2.2,1 Locomotive Parameters 

In Che methodology, a non-linear model is used to 
represent the locomotive. Figure 2 schemc^tlcally delineates 
the 15 degrees of freedom represented in the locomotive 
model. The model is designed to be as simple as possible, 
and yet to include all important suspension characteristics 
derived from the tests performed by Martin Marietta on sample 
locomotive trucks. 

Some important suspension features highlighted by the 
experimental test data were! 

- vertical friction damping on the truck frame ped- 
estals proportional to tractive effort, to wheel- 
rail steering forces, and to lateral thrust collar 
load on the roller bearings; 

- external hydraulic dampers having asymmetric exten- 
sion-compression characteristics; 

- centerplate rotational friction proportional to 
centerplate load; 

- secondary lateral friction damping proportional to 
tractive effort; 

- spring bottoming at all interfaces. 

Features not measured in the tests but included in the 
model are as follows: 
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FIGURE 2 - NONLINEAR MODEL DEGREES OF FREEDOM 
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- Flange contact:, 


- Non-linear wheel-rail creep forces dependent on 
inncantaneoua wheel loads, 

- Uneqvinl wheel diameters on different axles, 

- Unequal diameters of wheels mounted on the same 
axle (tape si 3 ;e mismatch). 

The complete list of parameters is included In the 
software description (Appendix 3). 

2 , 2.2 e r a t ipji u 1 P a r a mot or a 

Parameters which normally vary during operation in- 
clude speed, tractive effort, lateral components of coupler 
force, and wind loads. For convenience, nominal track 
curvature and superelevation are also included ns operational 
parameters . 


?. 3 Track (loometry Defects 

Four types of track geometry defects are modelled; 

1. Cross-level 

2. Lateral (or alignment of both rails) 

3. Vertical (or surface of both rails) 

4. Track gage variations. 

For each of the types above, defects can be specified 
in the form of transient pulses, as shown in Fig. 3 . Any 
number of those pulses can be specified along a section of 
track* The basic waveform is a sinusoid having an amplitude 
and a wavelength characteristic. Combined defects can be 
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FIGURE 3 - BASIC WAVEFORM OF TRACK GEOMETRY DEFECT FOR VERTICAL, 
LATERAL, CROSS-LEVEL, AND TRACK GAGE IRREGULARITIES. 
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relation BETWEEN DERAILMENT COEFFICIENT (WHEEL 
L/V RATIO) AND TIME DURATION REQUIREMENT FOR 
FLANGE CLIMBING. 

(Extract from "Outline of 250 Km/h Test run for the 
new Tokaido line" by Y. Ishihara. J.N.R. bulletin) 
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speciEied by superposition oC defects In any of the four 
types listed above. As well, other waveforms can be obtained 
by superposition of different wavelengths (i.e, Fourier 
components) . 

Although not specifically programmed in its present 
form, minor modifications would allow arbitrarily defined 
track defects. 

2 . 3 OUTPUT RESULTS 

Since the model computes time histories of all wheel- 
rail and suspension forces, there is complete flexibility 
in the selection of measures of safety. 

However, Judging from the technical literature, only 
three measures of safety have received enough attention both 
theoretically and experimentally to Justify inclusion in 
the methodology at this stage. These measures result from 
three specific derailment (or track damage) mechanisms as 
follows : 

2.3,1 Climb 

The primary indicator of wheel climbing is the ratio 
of lateral to vertical load on an individual wheel. This 
is the oldest known mechanism for derailment . (4 ) Under 
normal operating conditions, the potential for flange climb- 
ing is practically nil when the wheel L/V ratio is less than 
0.8. However, because the actual wheel climbing action 
depends on other factors, such as flange angle, wheel-rail 
coefficient of friction, and wheel-rail angle of attack, 
higher values than 0.8 may be sustained without derailment . (10) 


One aspect considered in the analysis is the fact 
that higher values of L/V can he sustained without derail- 
ment provided the time duration is short. The generally 
accepted figure for time duration is 50 msec, as shown in Fig» 

4. In the methodology this effect is implemented by the 
use of a simple "running average" over 50 msec. Since 
this time duration is adjustable in the model, values other 
than 50 msec, can be specified. 

2.5.2 lateral Tr ack P anel S hift 

This mechanism consists of a shift of the track 
structure (r-^ils and ties) laterally with respect to the 
ballast. Although not a derailment mechanism in Itself, 
it is included in the methodology as a mode of track failure. 
Experimental work in this area indicates that the tendency 
for track panel shift depends on the ratio of lateral to 
vertical axle load or net "wheelset L/V ratio”. (11) 

Critical values for wheelset l/V ratio depend on many 
track design details such as tie type, size, and spacing, 
as well as ballast type and shoulder width. Track panel 
strength is also quite sensitive to the degree of ballast 
consolidation. Typical critical values vary from 0.3 on 
newly worked track, to 0.7 or more on well compacted track. (3,11) 

2.3.3 Cage w idening by Rai l Rotati on (Rail Rollover) 

When the track gage becomes excessively large, (typ- 
ically 3| inches wider than the standard track gage), it is 
possible for a wheel to drop inside the rail onto the ties. 

Visual observations of the remaining trackage after derail- 
ments have sometimes revealed this condition. 
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Thcof etlenlly , two mechanisms can produce the lateral 
railhead motion required to obtain the wide 8»8®! first, 
lateral shift of the rail base and secondly roll of the 
complete rail section as a rigid assembly (neglecting small 
distortions of the rail section). 

Very little experimental data exists on the critical 
levels of load required to produce wide gage under typical 
"in-service" conditions. The problem is quite complex when 
considering the simultaneous action of many wheels of a 
truck under the effect of rail torsional stiffness. This 
is an area of current research by the AAR. (5,9) 

There is, however, one aspect that must be considered 
as a potential track damage mechanism. This is the total 
lateral to vertical load ratio exerted by a truck on one 
rail called "Truck Side L/V Ratio". Under high lateral loads, 
the resultant force vector on the rail can point outside 
the edge of the rail base. From this point on, rail roll- 
over is only prevented by the restraint from the spikes. 

Critical values for derailment depend on the "height 
to width ratio" of the rail section, as well as on the strength 
of the spikes or fastening system. In practice, an accept- 
able truck side L/V ratio would vary between 0.5 for an un- 
restrained rail to ( 0.5 + 20, 000/wheel load (lb) ) for 
wooden ties and new spikes(3). 
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l.lk SOFTWARE OVERVIEW 


The meuhodology uses a software package written In 
FORTRAN. The program is named PSE~2, for "jDynamics of 
^Jx-axle Locomotives", and was developed on a GDC Cyber 
174 computer. 

Many features of the program were designed to be user" 
oriented . Input parameter names and units were chosen to 
be those normally used by the railroad industry. The program 
operates interactively on a remote terminal with respect to 
input and output data. 

To facilitate the handling of the locomotive parameters, 
of which there are many, predefined or "canned" sets are used 
in the program and a facility is provided to alter any para-' 
meter under interactive control from the terminal. 

All input parameters are "echoed" for checking purposes, 
and all output formats are labelled with respect to date, 
time, run number Identification, and program revision number. 

Concerning software structure, the program is completely 
modular, containing over forty subroutines which provide 
flexibility with respect to program changes. 

Frogram running costs are reasonable, considering 
the complexity of the model. On the GDC Cyber 174 computer, 
execution requires a maximum of 11 seconds of CPU time for 
each second of real time simulation. Program execution 
requires 18K words of memory. Typical running costs on a 
commercial time sharing system vary between one and two dollars 
per second of real time simulation. 

More detailed software documentation is provided in 
Appendix III, 


CHAPTER 3 


APPLICATIONS. RESULTS. AND DISCUSSION 


3.1 INTRODUCTION 


The raethadology can play an important role as a pre- 
dictive technique, as well as a test support tool in many 
areas of rail safety. Typical applications are; 

- assessment of the importance of specific suspension 
design details; 

- comparison of different locomotive designs; 

- determination of approprlace maintenance standards 

on locomotive suspension elements (e»g. shock absorbers); 

- determination of acceptable track geometry defects 
and minimum track strength requirements; 

- Investigation of specific derailment mechanisms. 

In order to demonstrate how the methodology can be used, 
two distinct analyses were performed as follows. The first 
analysis consisted of an extonslvo parameter sensitivity study 
on a selection of twenty parameters relating to suspension 
design, operational conditions, and Crack geometry defects. 

The second analysis consisted of comparing the relative safety 
performance of three locomotive designs subjected to severe 
track geometry defects. Suspension parameters defining the 
three locomotives are based on Martin Marietta test data for 
the SDP-40F, the U-30G, and the E-8 locomotive trucks . (12 , 13 , 14) 
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3.2 PARAMETER SENSITIVITY ANALYSIS 


3.2.1 Do s cri p I i on 

In this study, a ref^erenee or "Base Case" la selected, 
fixing all required Input parameters. Evaluation of the safety 
performance is made by recording the peak or maximum value of 
each of the three measures of safety. (It is emphasized 
here that only one peak value is retained for each measure 
of safety without regard to the wheel or wheelset position 
at which it occurred.) Parameters are then changed one at 
a time, to determine their effect on the safety of operation. 

The reference case is chosen to be the SDF-40F locomotive 
equipped with "New HTC"* trucks, running on a 3-degree curve 
at a normal passenger speed of 65 mph. 

A single cycle of track geometry defect is used, contain- 
ing deviations in cross-level, alignment, and surface, and 
whose severity is set at the acceptable limit for class 4 

track on curves. The wavelength of the defect is chosen to 

n ) 

be 78 ft., similar to that of the "Perturbed Track Tests". 

Having defined all input data, a set of twenty parameters 
were selected for the sensitivity analysis. 

3.2.2 ResuJ t s 

The observations listed below represent the most 
significant findings from the sensitivity analysis (related 
to the new IITC truck design). However, most observations 
are of a sufficiently general nature that they probably apply 
to other locomotive designs. Details of procedures and numer- 
ical values are given in Appendix 1, 

^Modified HTC design (>Jan 1977).* Soft rubber bolster springs, 
lateral (secondary) shock absorber, steel pedestal liners (see 
Table 2.1). 
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a) The severity of the track geometry defect Is the most 
Important factor governing the peak values of the 
measures £ safety, Compared to those obtained for 
steady state values on a 3 degree curve (wheel L/V = >29, 
wheelset L/V “ .07, truck side L/V = .18), the presence 
of a severe geometry defect significantly increases 

the peak values. The safety measure which Increases 
roost drastically due to the presence of the defect 
is tlie wheelset L/V representing the potential for 
lateral track panel shift. 

b) Of the three components Included In the composite 
track defect used (i.e, surface, cross-level and align- 
ment) , the alignment component produces the largest 
peak values of the safety measures when used alone, 

c) For a given amplitude of defect, the wavelength has a 
large Influence. Of all the wavelengths studied (between 
39 ft and 118 ft), the shorter wavelengths give the 
largest peak values of safety measures, 

d) Wheel rail adhesion has a large influence on "wheel 
L/V" and "truck side L/V”, but little effect on "wheel- 
set L/V". 

e) Wheel size differentials have an effect on the safety 
measures. However, small mismatches In wheels of the 
same axle (i.e. from left to right) have a larger effect 
than differences between axles of a truck. 

f) Conicity of the wheel treads has a large Influence on 
the safety measures. Increased conicity (e.g. worn 
wheels) can reduce significantly the "wheel" and 
"truck side" L/V ratios. However, it was noticed that 
the peak value of wheelset L/V increases and begins to 
occur on the trailing axle of the trailing truck, 
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g) No individual auspenaion damping eiemant was found 
to have a major effect by itself, this is probably 
because there are several sources of damping on the 
new HTC truck, which all contribute to the overall 
system (e,g. pedestal friction, bolster friction, 
external vertical dampers, and lateral hydraulic 
dampers) . 
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3.3 COMPARATIVE ANALYSIS 


3.3.1 pegcTj P i 0 n 

The objective of this study was to demonstrate how the 
methodology can be used to compare different locomotive 
designs. Such comparisons can be used, for example, to 
highlight problem areas and possible solutions on existing 
locomotives, or to examine at the design stage the projected 
performance of new locomotives. 

In this analysis, three locomotives are compared, for 
which suspension characteristics are available from the Martin 
Marietta Corp. truck test data (12 , 13 , 14) . The three loco- 
motives are the SDP-40F equipped with New HTC (soft) trucks, 
the U-30C, and the E-8 locomotive. In order to provide 
realistic operational conditions, the following operational 
parameters are used for the simulated test. 

- 3 degree curve, 

- 6 inches superelevation, 

- speed range 40 - 65 mph^ and 

_ slnf.le track geometry defect at the 

allowable limit of FRA Class 4 track. 

Similarly to the sensitivity analysis, a combined 
defect Is used, having components In cross-level, alignment 
and surface. 
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In order to show possible dependence of the cowparative 

perforraanec on the type of track defect, two different 

wavelengths of track defect are used, The first 

n \ 

is set at 78 ft, sirailar to the PTT ' ^ procedure, and the 

» 

second is set at 39 ft, representing a standard rail length. 

Comparative performance of the three locomotives Is 
measured by the maximum values of each of the three measures 
of safety for the speed range considered, 

3.3.2 Results 

It must be emphasized at this point that the comparative 
analysis results only represent locomotive performance under 
a limited range of operating conditions. As well, the results 
only show maximum or peak values of the measures of safety 
and do not necessarily represent average or steady state values. 

Since the model has not been test verified, the absolute 
values of measures of safety may be questionable. however, the 
relative effects of parameter variations are considered valid. 

Details on procedures and numerical results are presented 
In Appendix 2. However, the following observations were made 
from the comparative analysis resul 1. 

a) Peak wheel l/V ratios (i.e. flange climbing) for the 
three locomotives are not substantially different 
from each other, and are generally lower than cur- 
rently accepted critical values for derailment. 

b) Peak wheel L/V ratios do not increase significantly 
as the speed is varied from 45 to 65 mph. 
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c) Significant difforences are found between the three 
locomotives when considering wheelset L/V ratio 
(i.ei track panel shift) resulting from the long 
wavelength defect considered (78 ft.)» when operat- 
ing at 65 mph (i.e. above the equilibrium speed). 

At this speed the U-30C locomotive showed a signifi- 
cant increase in wheelset L/V, when compared to the 
SDP-40F and E-8. However, simulations indicate that 
substantial improvements can be made to the U-30C by 
the addition of extra dampers in the secondary lateral 
S!>spension (i.e. between the truck and car body). 

d) Operating speed has a major effect on the peak 
wheelset L/V ratio towards the outside rail (i.e, 
track shift outwards on curve). Above the equilibrium 
speed, peak wheelset L/V ratios increase relatively 
sharply on the U-30C and the SDP-40F, but less on 

tihe E-8. 

e) Simulation results Indicate that the peak wheelset 

L/V ratio resulting from maximum allowable class 4 defects 
reach the critical values for unconsolidated track 
conditions . 
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3,4 COMMENTS ON YUE METHODOLOGY 


The purpose of this section is to discuss briefly 
various aspects considered In the methodology as they 
reflect real life conditions. As well, it is desired to 
highlight various modelling aspects regarding track and 
locomotive characteristics as well as derailment safety 
which are considered desirable as enhancements, This will 
help to provide a better understanding of the scope and present 
limitations of the approach, 

In the study of derailments, the investigation of 
transient response is more pertinent than that of steady 
state response because isolated large amplitude track defects 
are more representative of actual truck conditions which cause 
derailments. Although it is possible to represent measured 
track data, the use of idealized waveforms allows a better 
understanding of the plienomena Involved. Furthermore, no 
information is available on the types of track defects which 
actually led to the reported derailments , (2) 

The present model features most of the important 
locomotive suspension characteristics covered by the test 
data. There are a few items however for which there is no 
experimental data. Although not believed to be critical 
at this stage, it would be desirable to obtain experimentally 
validated data. These are: 

a) roller bearing lateral friction; 

b) measured ovortravel stop characteristics; 
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c) csCEccca 00 side-shearings on centreplace rotational 
8 tl 0 f ness I 

d) locomotive body torsional stl££neas; 

e) non-linear rubber characteristics; 

f) coupler stiCCneas and damping characteristics 
with respect to lateral and vertical motions. 

Considering the measures oC safety, the approach has 
been to compute peak values of the l/V ratios. TtiC user 
Chen has to make a Judgement on the actual potential for 
derailment or track damage. It would be desirable to build 
into the metho‘lology a set of safety criteria, indicating 
whether the level of safety is acceptable or not. To do 
this, a flange climbing model would have to bo incorporated 
in the model. To predict track damage, the algorithms would 
require changes because the critical values of many of the 
measures of safety are dependent upon the level of Che vert- 
ical wheel load. 



CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 

A methodology has been developed to compare the safety 
of operation of different locomotive designs. Its main 
elements are a method of characterization of track defects, 
a locomotive model with corresponding data on suspension 
characteristics, and a method of characterizing operational 
safety. The methodology has been used to study the effect 
of various locomotive parameters and operational conditions 
as well as to compare the operational safety of three different 
locomotives . 

Judging from the results to date, the methodology 
gives a faithful representation of locomotive behaviour. 

It has been possible to duplicate analytically many trends 
observed during previous experimental tests. It has also 
been possible to simulate little-understood behaviour observed 
in actual locomotives such as the occurrence of large dynamic 
forces on the trailing axle of the trailing truck. It Is 
evident that the approach taken is worth pursuing. 

The methodology can be used in other areas of rail 
safety not specifically investigated in this report. It can 
be used for example in refining information on acceptable 
track geometry defects, whether the defects are taken In- 
dividually or in combinations; it can also be used in the 
determination of safe operating speeds. It is believed that 
the methodology can play an important role in the investigation 
of derailments and in the study of specific derailment 
mechanisms. Another potential application is in the plan- 
ning of field experimentation, including the selection of test 
conditions, Che design of instrumentation procedures, and 
thedevelopmentofsafeCycriteria. 
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Although the present state of development of the 
methodology makes it a readily usable tool, many areas 
leave room for improvements, as discussed in the previous 
section. Briefly, improvements are possible in the track 
characterization, including the use of actual track defects, 
in the addition of other suspension features (e.g, side bear** 
ings), and in the definition and inclusion of safety criteria 

Although exhaustive checks have been conclucted to 
insure model validity, validation tests should be run to 
define simulation accuracy. data is recommended 

for use in this validation. Validation is important to 
ensure that significant factors have not been neglected 
in the design of the model and in the definition of the 
measures f safety. It is emphasized that a methodology 
validation shoiuld ideally extend into a range of operating 
conditions at or beyong limits of safety, 
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APPENDIX 1 


PARAMETER SENSITIVITY ANALYSIS 
1,1 SCOPE OF THE STUDY 

An analysis was made to determine the sensitivity of 
the three measures of performance (L/V ratios) to various 
locomotive and track defect parameters. Of the three loco- 
motives under study (Appendix 2), the SDP-40F with new HTC 
trucks was selected for sensitivity analysis since it is more 
common than either of the others. 

Except where otherwise noted, a standard track defect 
and set of operational conditions were used for the analysis. 
The defect, shown in Table 1.1, is a combination of crosslevel, 
lateral, and vertical misalignments having amplitudes equal 
to the maximum limit of FRA Class 4 track standards. This 
composite defect excites all modes of the dynamic system. 

The wavelength of the defect was chosen to be 78 feet, which 
is the same as that used in the Perturbed Track Test (PTT) . 

The stajidard operating condition was selected to be: 

i) 65 mph locomotive speed, 

ii) 3 degree curve with 6 inches superelevation, and 

iii) 9000 pounds of tractive effort. 

Statistically, it has been found that many derailments 
involving SDP-40F locomotives have occurred at medium speeds 
on 3 degree curves. (2) With 6 inches of superelevation, the 
maximum allowable vehicle speed on this curve would be 65 
Mph based on the conventional railway practice of a maximum 
of 3 inches cant deficiency. The 9000 lbs. of tractive effort 
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was reprcsencative of what might be required to pull a 
conventional passenger train at 65 mph through a 3 degree 
curve. 


The parameters selected for the analysis were those 
which could vary in actual service (for example: wheel/rall 
adhesion, defect amplitudes, flange clearance), those which 
might be affected by maintenance practice (for example: center 
axle diameter n^ismatch, external damper limiting forces), and 
those which could be changed with minor redesign to the equip- 
ment (for example: spring stiffnesses and travels). 

The ranges of the parameter variations were selected 
as those which could be found in railway service, or those 
which might bo achieved within the existing design. For 
example: flange/rail clearances could vary from 0.30" to 1,50", 
either because of deviations in track gauge which are within 
the acceptable FRA track specif Icatlons , or because of in- 
tentional gauge widening in curves. 

A total of 21 parameter sensitivities were evaluated 
which fall into one of the following six categories: 

wheel/ rail Interface 
wheel set geometry 
primary suspension 

- secondary suspension 
track defect geometry 
operating conditions 

The analysis was made by changing only one parameter, 
or only one set of parameters which in combination produce 
a single change to the System, at a time. For example: to 

study the influence of degree of curvature only, both the 
curvature and superelevation parameters were varied together, 
so as to have the same net lateral Eorce due to centrifugal 
acceleration acting on the locomotive. 
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1.2 


SENSITIVITY OF SELECTED PARAMETERS 


Tha following section pi^esents the results and a brief 
discussion of the parameter sensitivity analysis performed. 

It should be noted that the results pertain to the evaluation 
of safety measures (L/V ratios) and consequently, the trends 
indicated may not be exploitable in actual locomotive designs 
because of other considerations such as the vehicle ride quality, 
hunting stability, or wear. 

The results are given in graphs (Figures 1.1 to 1.21), 
with the three L/V ratios (individual wheel, wheelset, and 
truclc-side ratios) plotted versus the parameter being studied. 

The L/V ratios shown in figures were the maximum values which 
occurred in a given simulation run and are not specifically 
related to a particular wheel or wheelset. However, in general, 
it was found that the highest wheel L/V occurred on the outer 
wheel of the lead axle of the front truck; that the highest 
wheelset L/V was on the leading axle of the rear truck; and 
that the highest truck- side L/V was found on the outer side 
of the rear truck. 

Because only one parameter was varied at a time, the 
results sometimes showed only a small sensitivity to that 
parameter, whereas the sensitivity might be greater for a 
different set of conditions or for other values of the other 
parameters. (e.g. lateral hydraulic dampers in the secondary 
are more effective when the tractive effort is small as shown 
in Fig. 1. 22) . 
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Tables 1.2 and 1.3 are a qualitative sumtimry of the 
results of the parametric analysis. The sensitivity of each 
parameter is described as being either small, moderate, or 
large for each of the three L/V ratios. "Small" sensitivity 
is defined when the variation between the maximum and minimum 
L/V ratios is less than 0.05; "moderate" sensitivity is between 
0.05 and 0.10; "large" sensitivity is for variations greater 
than 0.10. 


1.3 


DISCUSSION ON PARAMETER SENSITIVITY 


1 • 3 . 1 Wheel/Rail Interface Parameters 


Three parameters in this group were studied: Flange/ 

rail stll£ness» flange clearancei wheel/rail adhesion. 

The flange stiffness which attenuates the lateral forces 
from the wheels to the ground, had only small effects on 
the I/V ratios. This is because the primary and secondary 
lateral suspensions of the locomotive are much softer than 
the range of rail stiffnesses considered, and hence the veh- 
icle stiffness has a predominant influence on the L/V ratios. 

Flange clearance influences the rolling radius difference 
between wheels and angle-o f-attack of the wheelset, and as 
a result, the creep forces. By increasing the flange clearance, 
the creep forces increase which tends to steer the wheelset 
to a more radial position and thus to decrease the L/V ratios. 

Wheel/rail adhesion is directly related to the creep 
forces, and provides a strong influence on the steering of 
the wheelset and hence the wheel L/V ratios. 

1.3.2 Wheelset Geometry Parameters 

Wheel conicity, center-axle mismatch, and wheel slde- 
to-aide mismatch were studied in this group of parameters. 

Both wheel conicity and wheel slde-to-slde mismatch have 
large influences on the L/V ratios, which is due to these 
parameters having a direct effect on the creep forces at 
the wheels. 
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Ccsncer-flxlii mlBmatch primarily changes the vertical load 
distribution within the truck, and has little effect on 
the steering moments of the truck* As a result, the L/V 
ratios do not vary significantly with center«axXe mismatch* 

1.3,3 Primary SuBpen&lon Parameters 

Five parameters were studied in this group: primary 

vertical travel and spring stiffness, pedestal friction, the 
limiting force of the external primary suspension dampers, 
and the lateral free play between the axle and journal boxes. 

It was expected that the h/V ratios would be affected by the 
first four of these parameters as a result of any vertical 
wheel unloading. Pot example, suppose a wheelset experiences 
a downwards vertical dip in the track, If Che wheelset is 
restrained from following the dip due to the primary 8uspenslon» 
particularly for the case of suspension bottoming, off-loading 
of the wheel/rail vertical force would occur. 

The sensitivity analysis for these parameters showed 
only small or moderate changes in the L/V ratios. This is 
because no substantial vertical off-loading occurs as a result 
of the track defect for the range of parameters selected. 

In addition, there is always some vertical damping present 
in the system (either due to pedestal friction or to external 
dampers, or both) which tends to control suspension-bottoming. 

Only small sensitivities of the L/V ratios were found 
for variations in the lateral free play between the axle and 
Journal boxes, This is because only a small amount of energy 
is removed from the system as a result of the friction between 
the axle and rollers in the Journal boxes and the lateral free 
play. 


1-6 



1 . 3 , 


Secondary Suspcnaxon Parameters 


The four parameters studied in the secondary suspension 
were lateral secondary travel and spring stiftness, the 
limiting force of the external secondary damper, and the 
bolster yaw friction factor, The first three parameters 
are related directly to the lateral dynamic forces which 
result from track defects, and showed moderate effects on 
the L/V ratios. These effects were noticeable at low values 
of the parameters, where suspension-bottoming and thus large 
lateral forces can occur. 

The bolster yaw friction factor principally Influences 
the truck swivel angle and hence the angles-of -attack of the 
wheelsets. For the range of friction factor studied, there 
were no appreciable changes in these angles and hence in the 
L/V ratios. 

1.3.5 Operational Parameters 

Two operational parameters were varied: degree of 

curvature and tractive effort, As discussed earlier, curva- 
ture and superelevation parameters were changed simultaneously 
such that there was no net change in the lateral force on the 
vehicle due to centrifugal forces. The Influence on the L/V 
ratios would then only be due to the degree of curvature. 

The results showed a small trend of increasing L/V ratios 
with degree of curvature; this is due to the increasing 
angles-of-attack of the wheelsets. However, because the 
curvatures were large, the trucks freely negotiated the curves 
without flanging of the trailing axles, and hence the angles- 
of-attack did not change significantly. 
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Tractive effort principally provides suspension damping 
due to friction in the pedestals and secondary suspension 
traction stops. The effects on the L/V ratios are small, 
as there are other sources of damping present in the syMtem 
which control the motions of the wheelsets and secondary 
suspension . 

1.3.6 Track Defect Geometry Parameters 

The amplitude, wavelength, number of track defects, and 
the type of track defects were studied in this group. 

f 

It was found that the L/V ratios are highly sensitive 

to variations in these parameters. This is .to be expected, 

\ 

since these parameters are directly related to the severity 
of the system forcing function. 
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I’RACK GKQMKTRY DEFECT CllARACTERISTlCS 


DRFECT No. 1 


COMPONENT 

TYPE 

AMPLITUDE * 
(in) 

WAVELENGTH 

(ft) 

1 

CROSS -LEV EL 

-1,25 

78 

2 

LATERAL 

1.5 

78 

3 

VERTICAL 

-1.375 

78 


TABLE l.l 


DEFINITION OF THE TRACK GEOMETRY DEFECT USED 
FOR TUE PARAMETER SENSITIVITY ANALYSIS. 


*Note: The following gives an equivalent description of 

the composite track geometry defect, in terms of 
the individual rails. 

Vertical, outer rail - 2 inches down, 

Vertical, inner rail ~ 3/4 inch down. 

Lateral, both rails - li inch to outside of curve. 



PARAMETER* 


SENSITIVITY // 




WHEEL 

L/V 

WHEELS ET 
L/V 

TRUCK-SIDE 

L/V 

1. 

Flange/rail 

stiffness 

Small 

Small 

Small 

2. 

Flange 

clearance 

Moderate 

Small 

Small 

3 , 

Wheel/rail 

adhesion 

Large 

Small 

Moderate 

4, 

Centre-^axle 

mismatch 

Moderate 

Small 

Small 

5. 

Wheel mismatch 
(right-to-left) 

Large 

Large 

Large 

6 . 

Wheel conicity 

Large 

Large 

Large 

7 . 

Primary vertical 
stiffness 

Moderate 

Small 

Small 

S . 

Travel-primary 
vertical suspension 

Small 

Small 

Small 

9. 

Lateral primary 
free play 

Moderate 

Small 

Small 

10. 

Pedestal friction 

Small 

Small 

Small 

11. 

Limiting force 

Small 

Small 

Small 


primary vertical 
dampers 


TABLE 1.2 - SUMMARY OF PARAMETER SENSITIVITY 
(page 1 of 2) 

* The parameter numbers correspond to Figures 1.1 through 1.21 
repsec tively . 

Smal3 A L/V <0.05 

Moderate ... A L/V » 0.05 - 0.10 
Large ...... A L/V >0,10 
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PARAMETER 


SENSITIVITY 




WHEEL 

WHEELS ET 

TRUCK-SIDE 



L/V 

L/V 

L/V 

12. 

Lateral secondary 
Stiffness 

Moderate 

Small 

Moderate 

13 . 

Travel-secondary 
lateral suspension 

Small 

Small 

Moderate 

14 . 

Limit force - 
external secondary 
lateral damper 

Small 

Moderate 

Small 

15. 

Bolster yaw 
friction factor 

Small 

Small 

Small 

16 . 

Tractive effort 

Small 

Small 

Small 

17 . 

Curvature 

Moderate 

Small 

Small 

18. 

Amplitude of 
track defects 

Large 

Large 

Large 

19. 

Wavelength of 
track defects 

Large 

Large 

Large 

20. 

Number of pulse 
of track defects 

Small 

Large 

Moderate 

21 . 

Type of defect 





: crosslevel 
: lateral 

Moderate 

Small 

Small 


alignment 
: vertical 

Large 

Large 

Large 


alignment 

Small 

Small 

Small 


TABLE 1.2 - SUMMARY OF PARAMETER SENSITIVITY (CONTINUED) 
(page 2 of 2) 
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APPENDIX 2 


COMPARATIVE ANALYSIS 


2 . 1 INTRODUCTION 


This section describes the selected input parameters, 
track geometry defects, and operating conditions used in 
the comparative analysis of the three following locomotives: 

1. SDP-40F locomotive with new HTC trucks 

2. U-30C locomotive 

3. E-8 locomotive 

The performance comparison is based on peak values of 
the three measures of safety, (i.e. wheel, wheelset, and 
truckside L/V ratios) when the locomotives are subjected to 
the same input functions. 

2 . 2 DEFINITION OF THE INPUT PARAMETERS 

2.2.1 Locomotive Parameters 

The input parameters for each locomotive are grouped 
as follows: 

1. Geometric Properties 

2. Mass Properties 

3. Suspension Properties 
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Table 2.1(a) gives a complete list of the Geometric 
Properties for all three locomotives. 

Table 2.1(b) defines the Mass Properties of the 
locomotives and trucks. 

The Input Parameters for the various suspension 
elements are defined in Table 2.1(c). Explanatory notes 
on some of the numerical values used are provided in Table 
2.1(d) . 

2.2.2 Track Geometry Defects 

A combination of crosslevel, lateral, and vertical 
misalignments having amplitudes equal to the ;maximum limits 
for FRA class 4 track safety standards are selected as track 
geometry defects. Two different cases are tested using the 
same defect amplitude, but with different wavelength. In the 
first case, the wavelength is chosen to be 78 ft (equal to the 
wavelength used in the Perturbed Track Test), and in the 
second case, 39 feet, which corresponds to the standard rail 
length. For both cases, each track defect is only used 
once (single pulse). Table 2.3 describes the details of each 
track geometry def^ict. 

2.2.3 Op era t Iona I Parameters 


The operational parameters are selected as follows: 

Track curvature is set at 3 degrees. This curvature 

has the advantage of satisfying the required range of track 

speeds for typical class 4 tracks. It also corresponds to 

the curvature over which many of the observed derailments 

( 2 ) 

occurred on six-axle locomotives. 
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Wifcli 6 inches of superelcvn tion , a mr .^mura of 3 Inches 
of superelevation deficiency, and using the p^ji'.ulure il- 
lustrated in Table 2,2, the minimum and maximum vehicle speeds 
were found to be 40 mph and 65 mph respectively. 

In summary, the selected operational parameters are; 

1. 3 degree curve, 6 inches superelevation. 

2. Speed range: 40 to 65 mph 

3. 9000 pounds tractive effort. 

The same parameters are used for both combined defects 
and for all three locomotives. 

2 . 3 COMPARATIVE ANALYSIS RESULTS 

Results for all simulations are shown in Table 2.4(a) 
and (b), indicating peak values of the measures of safety 
and the position on the locomotive where the peak values 
occurred. Results of peak values are also displayed graph- 
ically in Figures 2.1 to 2.6. 

2 • 3 • 1 Comparative Perfo rmance on the 78 ft Wavelength Defect 

Considering "Wheel L/V", indicating flange climbing 
potential (Fig. 2.1), only minor differences exist between 
the three different locomotives. As well, since the peak 
wheel L/y ratios are lower than currently accepted critical 
values for derailment (typically above 0.8), it appears that 
the probability of derailment on the selected defect due to 
flange climbing is small. 
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In Figure 2.2, relatively large differences are seen 
between the locomotives when considering wheelset 1/V, indicat- 
ing the potential for track panel shift. At 65 mph the U-30C 
produces peak values nearly twice as high as the otlier two. 

It can be noted as well that this peak value occurs on the trailing 
axle of the trailing truck. A similar response is found on 
the U-30C when considering truck side L/V (Fig. 2.3). 

2.3.2 Goiiiparative Performance on the 39 ft WavelonRt h hefect 

Peak wheel L/V ratios are not found to be significantly 
different on the three locomotives, the differences being in 
the 10 percent or less range (Fig. 2.4). The SDP-40P gave the 
lower values for all speeds. On wheelset L/V ratio (Fig. 2.5), 
the E-8 locomotive gives the higher Values, but the difference 
is more pronounced at lower speeds than at 65 mph. 

On truck side L/V ratio, the U-30C locomotive shows the 
best performance, contrary to what was shown previously for 
the 78 ft wavelength. 


2.4 


DISCUSSION 


The largest and probably most significant difference 
in behaviour observed between the three locomotives is the one 
resulting from the long wavelength combined defect at 65 mph. 

In this particular case, the U-30C locomotive is found to 
produce significantly higher peak values of wheelset L/V and 
truck side L/V, as compared to the other two locomotives. 

In order to provide an engineering explanation for this 
situation, the parameters of the U-30G are compared to those 
of the SDP-40F. The SDP-40F locomotive is chosen for reference 
because its weight and inertial characteristics are not sig- 
nificantly different to those of the U-30C. In this case, 
only the suspension characteristics can be responsible for 
such a large difference in behaviour. 

If one examines only those suspension parameters 
which show differences of 50 percent or more, three parameters 
are found as follows: 

1. Lateral secondary suspension dampers. 

2. Primary vertical stiffness. 

3. Pedestal friction. 

In order to quantify the effect of those three para- 
meters, additional runs were made with parameter modifications 
as shown in Table 2.5. 

It is clear from the numerical results of Table 2.5 
that the absence of lateral secondary damping in the U-30C 
locomotive is the explanation for its different behaviour as 
compared with the other two. 
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INPUT PARAMETERS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


DEFINITION 


GEOMETRIC PROPERTIES 

LOCOMOTIVE BODY 

COUPLER PIN/CG-LONG. 
SECONMRY/CG-VERT. 
PRIMARY/SECONDARY-VERT. 
COUPLER/RAIL-VERT. 

CENTER OF PRESSURE/RAIL-VERT. 
TRUCK CENTER/CG-LONG. 


UNITS SDP-40F U-30C 


TRUCK FRAMES 

HALF TRUCK WHEELBASE 
HALF LATERAL SPACING 
THICKNESS OF SHIMS 
IN THE VERTICAL 
PRIMARY SUSPENSION 

(POSITIVE UHEN 
SHIMS ADDED) 


LEFT AXLE 1 IN 

2 IN 

3 IN 

4 IN 

5 IN 

RIGHT I !n 

2 IN 

3 IN 

4 IN 

5 IN 

6 IN 



HHEaSET/TRACTION MOTOR ASSEMBLY 
NOMINAL TREAD RADIUS 
HHEEL TREAD CONICITY 
DIFFERENCE IN MEAN 
ROLLING RADIUS FROM 
NOMINAL 

(POSITIVE WHEN LARGER 
THAN NOMINAL) 

MISMATCH IN NHEa ROLLING 
RADIUS FROM SIDE TO SIDE 
(POSITIVE UHEN WHEEL ON 
LEFT RAIL IS LARGER THAN 
THAT ON RIGHT RAIL) 


HALF KINEMATIC GAUGE 
RAILHEAD CROWN RADIUS 


see notes cable 2.1(d) 


IN 

IN 

AXLE 1 IN 

2 IN 

3 IN 

4 IN 
B IN 
6 IN 

AXLE 1 IN 

2 IN 

3 IN 

4 IN 

5 IN 

6 IN 


20.000 

.050 

0.00 

0.00 

0.00 


20.000 

.050 

0.00 



TABLE 2.1 - INPUT PARAMETERS 


(a) geometric PROPERTIES 
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INPUT PARAMETERS 


DEFINITION 


MASS PROPERTIES 

LOCOMOTIVE BODY 
HEIGHT 

ROLL MOMENT OF INERTIA 
PITCH MOMENT OF INERTIA 
YAH MOMENT OF INERTIA 


TRUCK FRAMES 
HEIGHT 

ROLL MOMENT OF INERTIA 
PITCH MOMENT OF INERTIA 
YAH MOMENT OF INERTIA 


HHEaSET/TRACTION MOTOR ASSEMBLY 
HEIGHT 


ROLL MOMENT OF INERTIA 
ABOUT CG 


OFFSET OF CG FROM 
AXLE CENTERLINE 

(POSITIVE HHETJ MOTOR 
TRAILS AXLE) 

YAH MOMENT OF INERTIA 


AXLE 1 
Z 

3 

4 

5 

6 

AXLE 1 


AXLE 1 
2 

3 

4 

5 

6 

AXLE 1 . 


UNITS SDP-AO 


LB 

LB-IN-SECZ 

LB-IN-SEC2 

LB-IN-SECZ 


LB 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 


LB 

LB 

LB 

LB 

LB 

LB 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

IN 

IN 

IN 

IN 

IN 

IN 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 

LB-IN-SECZ 



161400. 

I6I400. 


13IZ4. 



131Z4. 

20000. 

20000 . 

ZOOOO. 


wiM'l 

/iliWtj 


10.0 

10.0 

10.0 
- 10.0 
-lO.O 
-lO.O 

16780.0 

16780.0 

16760.0 
16780.0 
16780.0 
16780.0 


U-30C 

E-8 

300700. 

ZIOBOO. 

1720000. 

1070000. 

39600000. 

Z5ZOOOOO. 

39600000. 

25200000. 

14510. 

15440. 

56000. 

55000. 

< 178000. 

170000. 

178000. 

1 

170000. 

11580. 

11580. 

11580. 

*1 3Z00. 

11580. 

11500. 

11580. 

11500. 

11560. 

3Z00. 

11580. 

11500. 

16750. , 

16000. 

18750. 1 

* 7000. 

16750. 

18000. 

1B7S0. 

16000. 

18750. 

* 7000. 

18750. 

18000. 

10.0 

10.0 

10.0 

0.0 

10.0 

-lO.O 

-10.0 

10.0 

-10.0 

0.0 

-10.0 

-10.0 

16500.0 

16500.0 

16500.0 

7000.0 

16500.0 

16500.0 

16500.0 

16500.0 

16500.0 

7000.0 

16500.0 

16500.0 


* see notes table 2.1(d) 


TABLE 2 . 1 (CONTINUED) 


(b) MASS PROPERTIES 
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OWGINAU PAGE jS 

OF POOR QUALITY input parameters 


DEFINITION 

UNITS 

SDP-40F 

U-30C 

E-8 

NO. 

SUSPENSION PfiOPESTIES 

SECONDARY 

LATERAL 

STIPfNESS 





1 

1 

1 

RUBBER STIFFNESS 
SPRING TRAVEL TO STOPIEACH SIDEI 
STIFFNESS BEYOND STOP 
DAHPINB ■ 

VISCOUS COEFFICIENT (RUBBER) 

FRICTION BOLSTER COEFFICIENT 

LVB COEFFICIENT FOR BOLSTER FRICTION 

LB/IN 

IN 

LB/IN 

ISSOO. 

1.Z50 

ISOOOO. 

12000. 

1.500 

150000. 

* 5000, 
2.2S6 
ISOOOO. 

s 

*71 

LB/(IN/SECI 

LB/(1N/SEC> 

* 136. 
.40 
4000. 

* 65. 

.45 

4000. 

‘ & 
4000. 

74 

YAH DAHPING 





ROTATIONAL FRICTION FACTOR 

IN-LB/LB 

1.10 

1.10 

l.IO* 

•75 

LVB COEFFICENT 

IN-SEC 

10000000. 

10000000. 

10000000. 

76 ! 

EXTERNAL DAMPERS ' 
VISCOUS COEFFICIENT 

LB/IIN/SECI 

lOZO. 

0. 

* 4000. 

” i 

LIMITING FORCE 

LB 

3600. 

0. 

*3000, 

76 ' 

! 

PRIMARY 

LATERAL 

STIFFNESS 






ROLLER BEARING FREE PLAY 

IN 

.31 

.31 

.31 

79 

RUBBER ELEMENTS TRAVa 

IN 

.2S 

.25 

.25 

60 

RUBBER PRaOAD 

LB 

1500. 

1500. 

1500. 

81 

STIFFNESS BEYOND STOP 
DAMPING 

FRICTION FOR ROLLER BEARING 

IN 

■ 

100000. 

100000. 

1 00000. 

*82 


.10 

.10 

.10 

83 

LVB COEFFICIENT FOR FRICTION 

LB/IIN/SEC) 

4000. 

4000. 

4000. 

84 

VERTICAL 

STIFFNESS 






AVAILABLE TRAVELICOMPRESSIONI 

IN 

Z.2S 

2.25 

2.25 

85 

AVAILABLE TRAVEL (EXTENSION) 

IN 

2.Z5 

2.25 

2.25 

BE 

SPRING STIFFNESS(PER JOURNALI 

LBMN 

5630. 

7920. 

*4900, 

87 

STIFFNESS BEYOND STOP(EITHER HAY) 

LB/IN 

200000. 

200000. 

200000. 

*BB 

PEDESTAL DAMPING 




FRICTION C0EFFICIENT(PEDESTAL FACE) 

— 

.35 

.16 

* 0.00 

89 

FRICTION C0EFFIC(ENT(SIDE LUG) 

— 

.35 

.16 

* 0.00 

90 

LVB COEFFICIENT FOR TOTAL FRICTION 

LB/(IN/SECI 

3000. 

3000. 

0. 

91 

EXTERNAL DAMPERS (OPTIONAL) 






CODE TO INDICATE THE LEFT AXLE 1 

— 

0. 

1. 

* 1. 

92 

PRESENCE OR ABSENCE OF 2 

— 

1 . 

0. 

1. 

93 

EXTERNAL DAMPERS 3 


0. 

1. 

1. 

94 

, < 


0. 

1. 

1. 

95 

(O'NO DAMPER) , . k 5 

— 

1. 

0. 

1. 

1. 

96 

(1=DAMPER) e 

— 

0. 

I. 

97 

RIGHT i 


0. 

1. 

1. 

98 

Z 


1. 

0. 

1 . 

99 

3 


0. 

1. 

1. 

100 

4 

— 

0. 

• • 

1. 

101 

5 


1. 

0. 

1. 

lOZ 

6 


0. 

t. 

1. 

103 

VISCOUS COEFFICIENT (COMPRESSION) 

LB/( IN/SEC) 

510. 

2000. 

3000. 

104 * 

VISCOUS COEFFICIENT(EXTENSION) 

LI/(IN/SECI 

510. 

3000. 

3000. 

105 

LIMITING FORGE(COMPRESSION) 

LB 

1800. 

1000. 

1500. 

106 

LIMITING FORCE(EXTENSION) , 


1600, 

1000. 

1500. 

107 

HHEEL/RAIL 

aANGE 

iN 





aASGEHAY CLEARANCE) EACH SIDE) 

.550 

.550 

.550 

108 

STIFFNESS AFTER FLANGE CONTACT 

LB/IN 

eoooo. 

^ flOOOO. 

80000. 

109 

TREAD 

CREEt COEFFICIENT ADJUSTMENT FACTOR 


.70 

.70 

.70 

no 

COULOMB FRICTION COEFFICIENT 

” . 

.30 

.30 

■ 30 

111 


* sue notes TABLE 2,1 (CONTINUED) 


(c) SUSPENSION PROPERTIES 
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table 2.1 (CONTINUED) 


(d) NOTES ON LOCOMOTIVE PARAMETERS 

Remarks 

Axles assumed equally spaced for U-30C 
and SDP-40F. 

22 Standard AAR wheel profile. 

46 Estimated, without traction motor. 

52-55 Estimated. 

69 Equivalent to theoretical stiffness of swing 

hangers (ref. 14, niain report). 

71 Estimated. 

72 Equivalent viscous coefficient at 2.5 Hz, based 

on test data at 0.25 Hz. For the E-8 locomotive 
see 77-78 below. 

77-78 Secondary lateral friction damping, independent 

of tractive effort as per test data. 

87 Equivalent to series-connected primary and 

secondary spring Stiffness. 

89-90 Test data indicates pedestal friction independent 

of tractive effort (see 92 to 107 below). 

92 to 107 Equivalent to measured test data. 
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Parameter No . 
7 


I 



ORIGINAL PAGE !3 
OF POOR QUALITY 


TABLE 2. 2 " SUPERELEVATION PARAMETERS 


Curve Paroroeter 8 : 3 degrees 

6” superelevation 


Formula ; 

Superelevation ? 
(Required ) 

D - 
V - 


0.0007 D 

degree of track curvature 
per 100 ft chord 

speed (m . p . h . ) 


Superelevation s 
(Def iciency) 


superelevation 
( Required ) 


superelevation 

(Actual) 


SPEED 

(mph) 

— M.| 

REQUIRED 

SUPERELEVATION 

(in) 

ACTUAL 

superelevation 

(in) 

SUPERELEVATION 

DEFICIENCY 

(in) 

40 

3. 36 

6 

-2.64 

45 

4.25 

6 

-1.75 

50 

5.25 

6 

-.7 5 

55 

6.35 

6 

. 35 

60 

7.56 

6 

1.56 

65 

8.87 

6 

2.87 
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TABLE 2.3 - TRACK GEOMETRY DEFECT CHARACTERISTICS 


DEFECT No. I 


COMPONENT 

TYPE 

AMPLITUDE* 

(In) 

WAVELENGTH 

(ft) 

1 

CROSS-LEVEL 

-1.25 

78 

2 

LATERAL 

1.5 

78 

3 

VERTICAL 

-1.375 

78 


DEFECT No. 2 


COMPONENT 

TYPE 

AMPLITUDE* 

(in) 

WAVELENGTH 

(ft) 

1 

CROSS-LEVEL 

-1.25 

39 

2 

LATERAL 

1.5 

39 

3 

VERTICAL 

-1.375 

39 


* Note: The following gives an equivalent description of 

the composite track geometry defect, In terms of 
the Individual rails. 

Vertical, outer rail - 2 Inches down. 

Vertical, Inner rail - 3/4 Inch down. 

Lateral, both rails - ij Inch to outside of curve, 
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tABLB 2.4 - COMPARATIVE ANALYSIS RESULTS 


(a) 78 ft Wavelength Geometry Defect 


LOCOMOTIVE 

SPEED 

(raph) 

PEAK VAL 

UES OF SAFETY IN 

DICATORS 



Wheel 

L/V 

Wheelset 

L/V 

Truck Side 
L/V 

SDP-40F 

40 

.51(L1,L4) 

,15(R2,R5) 

.18(Rl,Ll,R2) 


45 

.52(L1) 

.12(R2) 

.21(L1,L2) 


50 

,53(L1) 

.10(R2,L4) 

. 24(L2) 


55 

. 54<L1) 

,12(L4) 

. 28(Ll) 


60 

.55(L1) 

.16(L4) 

.30(L1) 


65 

.57 (LI) 

.19(L4) 

. 33(L2) 

U-30C 

40 

,5S(L1) 

.16(R2) 

. 22(L2) 


45 

.55 (LI) 

.14(R2) 

. 24(L2) 


50 

.55 (LI) 

.12(R2,R5) 

,27(L2) 


55 

. 54(L1) 

.11(L4) 

.28(L2) 


60 

.55(L1) 

.21(L4) 

.32(L2) 


65 

.55(L1,L4) 

.35(L6) 

.43(L2) 

E-8 

40 

.54(L1) 

.19(R5) 

. 22(L1) 


45 

.54(L1) 

. 15(R2) 

.24(L1) 


50 

.54(L1) 

.1S(R2) 

.26(L1) 


55 

.55(L1) 

. 12(L1,L4) 

. 28(L1) 


60 

.57(Ll) 

.15(L4) 

. 30(L1) 


65 

. 60 

.17(Ll,L4) 

. 32 (LI) 


Note: Designations in brackets denote position at which peak 

values occur on a right-hand curve. 


L, R ■ Left & Right sides. 

1 to 6 ■ wheel or wheelset, numbered from the lead axle. 
1,2 “ lead & trail, for truckside only. 


2-12 











table 2. a - COMPARATIVE ANALYSIS RESULTS (CONTINUED) 


(b) 39 ft Wavelength Geometry Defect 


LOCOMOTIVE 

SPEED 

(mph) 

PEAK VALUES OF SAFETY INDICATORS 

Wheel 

L/V 

- - 

Wheelset 

L/V 

Truck Side 
L/V 

sdp-aof 

40 

.74(L4) 

. 20(L1,R5) 

.27(L1) 


45 

.76(L4) 

,22(L1) 

.35(L2) 


50 

.78(L4) 

. 24 (L4) 

.39(L2) 


55 

. 80(L4) 

.30(l4) 

.40(L2) 


60 

.83(L4) 

. 34(L4) 

.40(L2) 


65 

.88(l4) 

. 3 9 (l4) 

.44(L2) 

U-30 

40 

. 75(L4) 

•22(R2) 

.23 (LI) 


45 

. 78(L4) 

. 21(R2) 

.25(L1,L2) 


50 

. 82(L4) 

. 23(L4) 

.33(L2) 


55 

.86(L4) 

. 29(L4) 

. 35(l2) 


60 

.87(L4) 

.34(L4) 

. 33(L2) 


65 

.87(L4) 

. 39(l4) 

.37(L2) 

E-8 

40 

. 76(L1 ,R2) 

.22 (LI 

. 27(L1) 


45 

.83(L1) 

. 29(L1) 

. 36(L1) 


50 

.84(L1) 

.32(L1) 

.42(L1) 


55 

.83 (LI) 

. 31(L1,L4) 

. 39(L1,L2) 


60 

.87 (LI) 

. 36(L1,L4) 

.39(L2) 


65 

.91(L1) 

.42(L1,L4) 

.41(L1,L2) 


Note; Designations in brackets denote poiiition at which peak 
values occur on a right-hand curve. 


L, R = Left & Right sides. 

1 to 6 = Wheel or wheelset, numbered from the lead axle. 
1,2 = lead & trail, for truckside only. 
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TABLE 2.5 - SENSITIVITY ANALYSIS ON THREE SELECTED 
PARAMETERS OF THE U-30C LOCOMOTIVE 


DESCRIPTION 

(U-30C Locomotive, 78 ft 
wavelength defect at 65 mph) 

MEASURES OF SAFETY 
(peak values) 

Wheel 

L/V 

Wheelset 

L/V 

Truck Side 
L/V 

NOMINAL CASE 

.55 

.35 

.43 

MODIFIED LATERAL 
DAMPING (1) only 

.58 

.27 

.35 

MODIFIED VERTICAL 
STIFFNESS (2)only 

.54 

,34 

.42 

MODIFIED PEDESTAL 
FRICTION (3) only 

.55 

.36 

.42 

(1) Addition of hydraulic da 
suspension, having 3600 

mpers in the lateral see 
lb limiting force. 

:ondar y 


(2) Primary suspension vertical stiffness reduced from 7920 
to 5630 Ib/ln per J</urnal. 


(3) Coulomb friction coefficient Increased from 0.16 to 0.33. 
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FIGURE 2.2 - PEAK WHEELSET L/V RATIO - 78 FT WAVELENGTH COMBINED DEFECTS 
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APPENDIX 3 


GENERAL SOFTWAkE DESCRIPTION 


This section describes the software used to implement 
the locomotive model. Presently, there are two programs, 
the first of which accepts the parameters defining a 
particular run with a particular locomotive, and performs 
the integration which constitutes the run. This program 
outputs most of its results to a disc file, which is used 
by the second program to plot selected variables versus 
time. This chapter deals only with the first or main 
program as it contains the model itself. 

The specification of input parameters is effected 
through two mechanisms. In the case of parameters 
definirg the three locomotives under study, default sets 
of parameters are specified within the program itself. 

All other input, including the specification of which set 
to use and the definition of track defects and of operational 
parameters (e.g. locomotive speed, nominal track curvature, 
etc.) is made by the operator via the terminal. 

Output from the program has three destinations. 

Some output goes to the terminal, especially as required 
for Interaction. A second group of data, including lists 
of all parameters specifying the particular run and the 
summary output of peak L/V ratios, is output to a disc 
file for subsequent (optional) listing, either to the terminal 
or to a high-speed printer. This file (file 116) may be 
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saved for future reference, Moat of the output from the 
program, however, is written to a second disc file (file //8) 
for further processing off-line. This information consists 
of blocks of data defining the state of the locomotive at 
fixed Intervals during the run. Presently, these data are 
used off-line to generate the plots of selected variables 
versus time. 

A heading identifying the run number and date and time 
of execution is written to the terminal and to both disc 
files. This heading ensures positive identification of data 
after a run. 

Except for a flowchart of the main features of the 
program, the descriptions of the software in sections 3.1 and 
3.2 are expressed In terms of the program structure. Anyone 
with a rudimentary appreciation of computer programs should 
have no trouble following the description of the structure 
and of the individual sof tware modules . 

3.1 OVERALL PROGRAM OPERATION AND STRUCTURE 


A flowchart showing the major operations of the program 
is given in Figure 3.0. 

At the beginning of a rurt the operator defines parameters 
relating to the vehicle, tr.'sck, and Integration which specify 
the run. Subsequently, other required parameters are calculated 
and the operator specifies the initial conditions for the 
locomotive. At this point, all the input parameters are 
reported and integration can begin. 

Each step of the integration entails the calculation 
of the items in the fourth box in Figure 3.0, the peak L/V 
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FIGURE 3.0 - COMPUTER CODE FLOW CHAAT 
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ratios (the measures of performance) are updated, and one 
step of Euler Integration is performod, Output of the 
variables describing the state of the locomotive is also 
performed, although usually less frequently than each 
integration step (this frequency is one of the operator- 
specified parameters). The above process is repeated until 
the run is complete, at which time a summary of the peak 
L/V ratios is output. 

These operations are described in more detail below. 

The balance of this section describes the structure 
of the major elements of the program which simulates the 
dynamic response of a locomotive. We can group the major 
tasks into five categories as follows: 

1. Specification of Input Parameters. 

2. Initialization prior to Integration. 

3. Reporting of Input Parameters. 

4. Integration and Output during the 
Integration Process. 


5. Summary reporting at the end of a run. 


These five basic categories are dealt with in detail 
in subsections 3.1.1 through 3.1.5 below. 

Figure 3.1 shows the structure (or hierarchy) of the 
major elements of the program. This figure Indicates that 
the main program (SMAIN) has control over twelve major 
elements which are grouped into the five categories mentioned 
above and labelled at the top of each column. In two cases, 
a second level in the hierarchy is shown; specifically, 
the routine SINVP has control over SPLIST and the integration 
algorithm (SEULER) has control over the routine 
responsible for providing the Instantaneous derivatives of 
the state variables (SCALGD) . 

Detailed descriptions of the functions of the individual 
routines in the five basic categories are contained under 
3.1.1 through 3.1.5. 

3.1.1 Specification of Input Parameters 

The first function of the main program is to allow the 
operator to specify the parameters required for a specific 
run. These parameters fall into three categories the 
input of which are handled by the following three routines; 

1) SINVP allows the operator to specify those parameters 

which describe the dimensional, inertial, and suspension 
characteristics of the locomotive to be modelled. 

This routine allows the operator to choose one of three 
sets of default parameters (over 100 parameters are 
involved) , and also allows him to specify individual 
parameters at his discretion. 

SINVp uses the routine SPLIST to list the parameters 

on the terminal if this is requested by the operator. 
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FIGURE 3.1 - STRUCTURE 





2) SXNTP allows the operator to specify the track 
defect definitions. Presently, the program allv>ws up 
to a total of 25 separate defects to be specified for 
any single run. Each defect is defined by its type 
(e.g. vertical, cross-level, etc. ) , its class (e.g. 
versed sine, ramp, step-, etc.), the position of the 
beginning of the defect (or defects) on track, and 
the wavelength, amplitude, and number of defects to 
be repeated. 

3) The third input routine, SINOP, allows the specification 
of operational parameters including: 

- Locomotive speed 

. - Nominal track curvature 

- Nominal track superelevation 
Total time for integration 
Time Interval for Integration 

Time Interval for output during integration 

- Locomotive tractive effort. 

The main program's only function, for this category, 
is to call the above three routines. 

3.1.2 Initialization Prior to Integration 

This group (the second column in Figure 3.1 ) provides 
for the initialization of parameters and allows the specifica- 
tion of the initial state variables. 

1) The routine SCALCDC calculates a number of secondary 
parameters which are derived from the parameters 
described in section 3.1.1. For example, the loco- 
motive velocity (specified in units of miles per hour) 
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is converted to units of feet per second and this 
converted value is stored for later computation, 

This routine is used to avoid the recalculation of such 
secondary parameters at each iteration. 

2) The other routine in this category is SXNIC, which 

is responsible for the balance of the initialization 
required prior to integration, The variables initialized 
fall into the following categories: 

a) State Variables 

The thirty state variables are all Initialized 
to zero and the operator is given the opportunity 
to modify as many as he wishes, 

b) Variables required for the first Step of Integration 

Some of the variables used in the calculations are. 
not calculated prior to their use within one 
integration step, and normally the values from the 
previous step are used. These variables must be 
initialized prior to the first step of integration. 
Specifically, these are: 

" the vertical wheel loads 

the longitudinal creep forces at the wheel treads 
- the primary forces at each axle box in the 
lateral and longitudinal directions 
a variable describing the rotational velocity 
of each of the locomotive axles. 


C'\-^ 
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c) Time Variables 

The time for Integration is initialized to zero 
and the print time is also Initialized. 

d) Sp eclal Error Flags 

Flvfa flags are provided in the software to indicate 
error conditions, and all are initialized to zero. 

At this time, only two of these are in use, and 
they flag Instances of wheel unloading and primary 
vertical unloading. 

e) Variables used to provide the Measures of Evaluation 

There are numerous variables which are used 
internally to the program in order to calculate the 
filtered L/V ratios used for evaluation of safety. 
These are Initialized here. 

The function of the main program for this category is 
to call the SCALCD and SINIG and to request the run number 
from the operator. 

3.1.3 Reporting of Input Parameters 

Prior to the beginning of the Integration, the main 
program reports the input parameters as part of the standard 
output. The current data and time are requested from the 
operating system, and a heading also containing run number 
and locomotive type code is output. This is followed by ^ 
the calling of the following three routines which list their 
output to a file for later use. 
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1) The routine SPLIST Is used to output all of the 
locomotive parameters set up through the use of 
SINVP. (Note that SPLIST appears twice in 
Figure (1); it is used both by SINVP as described 
under 3.1.1, and also by the main program). 

2) The routine STLIST is used to keep a record of all 
the track defects specified through the use of SINTP. 

3) SOLIST outputs the operational parameters: locomotive 

speed, curvature, nominal supisrelevatlon, length of 
track. Integration time step and printing time step, 
and tractive effort. 

The calling of these three routines concludes the 

reporting of input parameters. 

3 • 1 • 4 Integration and Associated Output 


Having completed the functions described above, all 
parameters are initialized and integration may begin. The 
main program is responsible for calling the routines in the 
fourth column of Figure 3 ,x to perform each integration 
step (SEULER) and to provide the output of variables during 
integration (SOUTPT and SPLOT) . 

SMAIN initially calculates the number of integration 
steps required, and then a loop is executed for each step. 
Within the loop, a check is made to see if the incremental 
time for printing has elapsed; if So, any requisite output 
is performed through SPLOT and/or SOUTPT (as of September 
17, 1980, SOUTPT is not in use). Also within the loop 
a call is made to the integrator (SEULER) , which performs one 
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step of integration and updates all the state variables 
and the time. SEULER uses SCALED to provide the 
Instantaneous derivatives of the state variables at 
every step. 

- Output routines: 

1) SOUTPT 

This routine is used to list various state variables 
and forces to a specified output device. This routine 
was originally used to print results to the terminal 
during integration as an aid in debugging; it is not 
presently in use and only SPLOT is used for output 
during integration. 

2) SPLOT 

This routine collects 260 variables which describe 
the state of the locomotive and outputs these to a 
disc file for later off-line processing. This 
processing presently consists of plotting selected 
variables versus time. The data saved on disc by 
SPLOT Includes most of the variable data in the 
program. 

- Integrator 

3) SEULER 

This routine calls the routine responsible for the 
provision of current values of the state variables 
(SCALED), performs one Euler integration step on 
the thirty state variables, and increments the 
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integration time before returning control to the 
main program. It la called at every integration 
step . 


The routine SCALCD la treated separately under 
section 3.2 below* 

3.1.5 Summary Reporting at the End of a Run 

At the end of a run,, when integration has proceeded 
for the specified time, some of the data is listed to the 
file used for printing which includes the error flags 
and the peak values and times of the evaluation measures 
(the L/V ratios). This having been effected by SUMOUT, 
the run is concluded, 


3 . 2 DERIVATI.VF.S OF STATE VARIABLES 


The routine SCALCD is responsible for the calculation 
of the derivatives of all the state variables. This 
routine makes use of several other routines as shown in 
Figure 3.2, which fall into five major categories as 
follows: 


1. Definition of instantaneous track defect inputs to 
the individual axlas. 

2. Calculation of the relative displacements (or rotations) 
and velocities at interconnections. 

3. Calculation of the relative forces (or moments) at 
the interconnections, 
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4 . 


Calculotlon of the absolute accelerations for 
each degree of freedom. 


ORIGINAL PAGE IS 

OF POOR QUALI1Y 


5. Calculation of measures of performance . 

These five categories are covered in subsections 

3.2.1 through 3.2.5 below, and correspond to the five 
columns in I- igure 3.2. 

3.2.1 Instantaneous Track Defects 

>1 . ■ * - m I ■ ■ 111 l — l ill .1 I ■III. 

Given the instantaneous position of the locomotive 
on track, SUPIN determines the displacement, velocity, 
and acceleration terms corresponding to all track defects 
under each of the six axles. There is provision for defects 
in cross-level, vertical alignment, lateral alignment, 
and gage. 

3.2.2 Relative Displacements and Velocities 

The routines in the table below are used to calculate 
the Instantaneous relative displacements (or rotations) 
and velocities at the various interconnections in the model, 
given the state variables (absolute displacements and 
velocities), and the Instantaneous track defects. 


ROUTINE NAME 

INTERCONNECTION 

direction 

SYW 

WKEEL/RAIL 

LATERAL 

SZP 

PRIMARY 

VERTICAL 

STHETW 

WHEEL /RAIL 

YAW 

SYP 

primary 

LATERAL 

STHETS 

SECONDARY 

YAW 

SYS 

SECONDARY 

LATERAL 
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The relntive dlspXacemencs and velocities are 
calculated once per step and are used to determine the 
Interconnection forces and moments. 

3,2.3 Forces and Moments at Interconnections 

The ten routines in this group calculate the forces 
and moments at the suspension interconnections. Some of 
these forces and moments are calculated from relative 
displacements (or rotations) and relative velocities at 
the interconnections and the constitutive relationships 
describing the suspension characteristics (for instance, 
vertical primary loads are calculated based on spring 
deflections and friction). Others are calculated on the 
basis of force and moment equilibrium (for instance, 
secondary pitch moments are calculated based on the 
difference between the primary vertical loads for the end 
axles of each truck). 

The following table indicates which routines calculate 
which forces and moments. 
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ROUTINE NAME 

INTERCONNECTION 

FORCE/MOMENT 

1) 

SFPX 

, ,1 . I 1 .. . M, J___. 

primary 

LONGITUDINAL FORCES 

2) 

SFPY 

primary 

LATERAL FORCES 

3) 

SFPZ 

primary and 

SECONDARY 

VERTICAL FORCES 

Vertical forces 

4) 

SFSY 

SECONDARY 

LATERAL FORCES 

5) 

SMSX 

secondary 

ROLL MOMENTS 

6) 

SMSY 

SECONDARY 

PITCH MOMENTS 

7) 

SMSZ 

SECONDARY 

YAW MOMENTS 

8) 

SFWZ 

WHEEL/RAIL 

VERTICAL WHEEL LOADS 

9) 

SCREEP 

WHEEL/RAIL 

LONGITUDINAL TREAD 
CREEP FORCES AND 
LATERAL TREAD CREEP 
FORCES 

xo) 

SFF 

WHEEL /RAIL 

LATERAL FLANGE FORCES 


1) SFPX 


This routine calculates the longitudinal pedestal forces 
which, because of pedestal friction, provide damping 
in the lateral and vertical primary suspension. 

Because there Is no primary longitudinal suspension 
in the model, these forces are calculated based on the 
following equilibrium conditions : 

a) The sum of the two primary longitudinal forces on 
any axle equals the nominal tractive effort for 
that axle, and 
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b) The difference between the primary longitudinal 
forces on an axle is due to a moment resulting 
from the inertial forces and 
wheel/rail forces acting on the wheelset. 

2) SFPY 

This routine calculates the lateral primary forces 
at each axle box due to free play, stiffness and 
damping between the axle and axle box, and pedestal 
friction. 

3) SFPZ 

This routine calculates the primary vertical forces 
at each axle box based on the nonlinear vertical stiffness 
characteristic Including steps, friction, and optional 
external dampers. The secondary vertical force acting 
on the vehicle from each truck is also calculated here, 

4) SFSY 

&FSY calculates the lateral secondary forces including 
a nonlinear stiffness term with stops and damping terms 
due to the rubber suspension, friction, and any external 
damper. 

5) SMSX 

This routine calculates the roll moment transmitted 
to the car body by the truck frame due to primary loads 
in the vertical and lateral directions. 
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6 ) 


SMSY 


TI 1 I 9 routine calculates the pitch moment 
transmitted to the car body by the truck frame 
due to primary vertical loads. 

7) SMSZ 

This routine calculates the yaw moment at the 
center plate due to friction* 

8 ) SFWZ 

This routine calculates the vertical wheel loads 
at each wheel from the primary vertical loads* the 
acceleration of the wheelset vertically (Including 
gravity), the rCill Inertia and acceleration of the 
wheelset and the lateral tread and flange forces. 

The model assumes no wheel-lift, and an error Is flagged 
if wheel-lift occurs. 

9) SCREEP 

SCREEP calculates the longitudinal and lateral 
forces between the wheel tread and the rail at each 
wheel based on the following: vertical wheel loads, 

the relative lateral displacement of the wheel and 
the rail, the lateral velocity of the wheelset, and a 
variable describing the slip velocity of each wheelset* 
This slip parameter Is also updated by SCREEP. 

10) SEF 

This routine calculates the lateral flange 
•force at each wheel, 
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3*2.4 Acceleratlona £ov each PcKree of Freedom 


These four routines contain the equations of notion 
of the system which are of the form 


a s F/m 
or 

a s M/I 

where a 5 acceleration, F * force, m 5 - mass, or 
a = angular acceleration, M ;s moment, 

I s moment of inertie* 

The routines are as follows; 


1) SVCLE 

This routine calculates the accelerations for the 
five locomotive body degrees of freedom. 


2) SDDP 

This routine calculates the lateral accelerations 
of the truck frames. 


3) SDDTHET 

This routine calculates the yaw accelerations 
of the truck frame and wheelset assemblies. 


4) SDDY 

This routine calculates the lateral accelerations of the 
six wheelsets. 
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The calculation of, theae accelerationB allows SGALCD 
to provide all the sote variable derivatives. For 
diaplaceinenta , the derlvativea arc aimply the corresponding 
velocity state variables; for velocities, the derivatives 
are the accelerations calculated above. 

3,2.5 Calculation of Measures of Performance 

The four routines which make up this group arc used 
to generate the L/V ratios which are used as measures 
of pcrforroance. The routines define Instantaneous lateral 
and vertical forces, generate the various instantaneous 
L/V ratios, filter this data, and preserve the values and 
times of peak values. This Information is presented in 
the summary output (section 3,1.5). 
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3.3 INPUT PROCEDURES AND OUTPUT FORMATS 


This section Is divided Into two parts. The first 
of these describes the options available to the operator 
for definition and modification of input data. The second 
includes sample output displays and describes their formats. 



3.3.1 Input Procedures 


As described briefly In section 3.1| there are 
four different types of Interactive Input between the 
operator and the program. 

These four shall be dealt with in this section, in 
the order in which they occur during program operation: 

1. Vehicle Parameters 

2. Track Defect Definition 

3 . Operational Parameters 

4. State Variable Initial Conditions 

1) Vehicle Parameters 

Vehicle parameters are defined by the operator from the 
routine SINVP. He is Initially requested to enter the number 
of a default set of parameters; presently there are three 
sets corresponding to an SDP-40F with the new HTC trucks, the 
U-30C locomotive, and the E-8 locomotive in that order. 


Having selected a default set, the operator has the 
option to list the parameters, change an individual parameter 
or parameters, or exit SIHVP. If he opts to list the parameters, 
he is subsequently again given the option to list, change, 
or exit. If he opts to change parameters, he identifies the 
parameter of group of parameters by number (see section 
3.3. 2.1) and specifies the desired parameter value. 

This approach to vehicle parameters relieves the 
operator from having to specify each of 111 parameters while 
retaining complete flexibility, 
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2) Track BeCect Definition 


Any combination of track defects may be Introduced 
Into the piece of track to be traversed by the locomotive 
through the specification of six data for each: 

a) Type 

The defect type indicates whether a specified 
defect is a defect in crosslevel (type 1), 
lateral alignment (type 2), vertical alignment 
(type 3), or gage (type 4). 

b) Class* 

The class of the defect indicates its shape. 
Presently the classes available are: versine 

(class 1), ramp (class 2), sinusoid (class 3), 
and step (class 4). 

c) Starting Position 

The beginning of the defect is specified in terms 
of feet along the track profile from the initial 
position of the lead axle. 

d) Amplitude 

This value is the total amplitude of the defect 
in inches. The sign convention is as follows; 




Type of defect 

Positive amplitude 


1. 

Crosslevel 

left rail elevated 

»• 

2. 

Lateral alignment 

track deflected to the left 


3. 

Vertical alignment 

track deflected upwards 

tt 

4. 

Gage 

wide gage 


*Note that this is not the FRA track safety 
standard classification. 
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ORIGINAL FAGE 13 
OF POOR quality 

For crosslevel defects » the crossleyel angle 
of the track Is calculated assuioing a nominal 
track gauge of 59i**. 

e) Wavelength 

This value is thd total length, in feet, of 
a single defect. It is npt meaningful for step 
defects . 

f) Number of defects 

The defect specified may be repeated, and this 
value Indicates the number of consecutive 
Identical defects. 

For example, thn specification; 

Type - 2 

Class - 1 

Start - 1.00 

Amplltu<^^ ^ P-5 
Wavelength - 39 
Number - 4, 

defines a grpup pf ^pur consecutive 39-foot verslne 
defects in lateral track alignment where each 
defect has an amplitude of i inch and the first 
defect starts 100^ from the starting Ip^catlon 
pf the lead aylp. 

Presently, the opeFetop may specify from zero to 
25 such defects interactively for each run, and 
defects may be superimposed even if they are of 
the same type- 
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3) 


0 p e r a C i 0 n a 1 Parameters 


The following parameters are specified interactively 
by the operator; 

- Locomotive speed - m.p.h. 

- Track curvature - degrees per 100 ft chord 

- Track superelevation - inches in 
Tractive effort - lb 

- Lateral wind load lb 

- Lateral component of lead coupler force - lb 
Lateral component of trail coupler force - lb 

- Total simulation time - seconds 

■> Integration step size - seconds 

Time interval for printing - seconds 

The last three are seconds of simulated timCi not 
seconds of computer time. 

4) State Variable initial Conditions 

Table 3^1 lists the thirty state variables 
corresponding to the fifteen degrees of freedom. As a 
rule, these are initialized to zero at the beginning of 
any run» but the operator may interactively redefine 
any of them. 
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If 

SYMBOL 

DESCRIPTION 

POSITIVE 

DIRECTION 

UNITS 

1 

^1 

1st Wlieelset lateral displacement 

left 

in. 

1 

^2 

2nd Wheelset lateral dieplacement 

left 

in. 

3 

^3 

3rd Wheelset lateral displacement 

left 

in. 

4 

>^4 

4th Wheelset lateral displacement 

left 

in. 

5 

^5 

5th Wheelset lateral displacement 

left 

in. 

6 

^6 

6th Wheelset lateral displacement 

left 

in. 

7 


Lead truck lateral displacement 

left 

in. 

8 


Trail truck lateral displacement 

left 

in. 

9 


Lead truck yaw rotation 

left 

radians 

10 


Trail truck yaw rotation 

left 

radians 

11 


Locomotive body lateral displacement 

left 

in. 

12 

Z 

Locomotive body vertical displacement 

up 

in. 

13 

a 

Locomotive body roll rotation 

right 

radians 

14 

B 

Locomotive body pitch rotation 

forward 

radians 

15 

T 

Locomotive body yaw rotation 

left 

radians 

16 

B 

1st Wheelset lateral velocity 

left 

in/sec. 

17 

B 

2nd Wheelset lateral velocity 

left 

in/sec. 

18 


3rd Wheelset lateral velocity 

left 

in/sec. 

19 

B 

4th Wheelset lateral velocity 

left 

in/sec. 

20 

H 

5th Wheelset lateral velocity 

left 

in/ sec . 

21 

^6 

6th Wheelset lateral velocity 

left 

in/sec. 

22 

Pi 

Lead truck lateral velocity 

left 

in/sec. 

23 

mm 

Trail truck lateral velocity 

left 

in/sec. 

24 

H 

Lead tiruck yaw Velocity 

left 

radians /sec. 

25 

®2 

Trail truck yaw velocity 

left 

radians/sec. 

26 

• 

Y 

Locomotive body lateral velocity 

left 

in/sec. 

27 

Z 

Locomotive body vertical velocity 

up 

in/sec. 

28 

* 

a 

Locomotive body toll yelocity 

right 

radians/ sec. 

29 

6 

Locomotive body pitch velocity 

fotvard 

radians/sec. 

30 

• 

Y 

Locomotive body yaw velocity 

left 

radians /sec. 


TABLE 3*1 - STATE variables 
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3.3.2 Output Formats 


As described briefly in section 3.1, there are 
five types of output from the program to the terminal 
and to the disc files. 

These five are presented in exhibits with accompanying 
discussion in the following subsections! 

1. Vehicle Parameter Output 

2. Track Defects Output 

3. Operational Parameter Output 

4. Summary of Results 

5. Output during Integration 

1) Vehicle Parameter Output 

Exhibit 3. 3. 2.1 shows the format of the heading 
output by the main program for run identification in the 
first three lines of the exhibit, and the balance of the 
exhibit shows the output of the vehicle parameters. 

This is identical to the interactive output available 
at the terminal during specification of the vehicle 
parameters, and the indices on the right are used for 
parameter modification. The list is divided into geometric 
parameters i inertial parameters, and suspension parameters. 
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ORlQIftAL rA\a,i lU 

OF POOR QUALITY 

PPOGBAN DSC-Z (REVISION OATEI AS OF SEPT.29TH, 1900) 
RUN I : 26Z OF aO/09/aO. 13.18.07. 

LOCOHOTIVE PARAHETERS SET NUNBER.* 1 


LOCOHOTIVE PARAHETERS SET NUNBERI 1 


DESCRIPTION 


UNITS 

VALUE 

INDEX 

CEOHETRIC PROPERTIES 
LOCOHOTIVE BOUT 

COUPLER PIN/CG-LONB. 


IN 

404.000 

I 

SECONDARY/CG-VERT. 


IN 

S0.200 

2 

PRIHARY/SECONDARY-VERT. 


IN 

7.500 

3 

COUPLER/RAIL-VERT. 


IN 

34.500 

4 

CENTER OF PRESSURE/RAIL-VERT. 

IN 

120.000 

5 

TRUCK CENTER/CG-LONG. 


IN 

276.000 

6 

TRUCK FRAHES 

HALF TRUCK HHEELflASE 


IN 

B2.200 

7 

HALF LATERAL SPACING 


IN 

39.500 

8 

THICKNESS OF SHIHS 

LEFT AXLE I 

IN 

0.000 

9 

IN THE VERTICAL 

2 

IN 

0.000 

10 

PRIHARY SUSPENSION 

3 

IN 

0.000 

11 

(POSITIVE NHEN 

4 

IN 

0.000 

12 

5 

IN 

0.000 

13 

SHINS ADDED) 

6 

In 

0.000 

14 


RIGHT 1 

IN 

0.000 

15 


2 

IN 

0.000 

16 



IN 

0.000 

17 


4 

IN 

0.000 

18 


5 

IN 

0.000 

19 


6 

IN 

0.000 

20 

HHEELSET/TRACTION HOTOR ASSEHBLY 

NOHINAL TREAD RADIUS 


IN 

20.000 

21 

HHEEL TREAD CONICITT 


IN 

.050 

22 

DIFFERENCE IN HEAN 

AXLE 1 

IN 

0.00 

23 

ROLLING RADIUS FRON 

2 

IN 

0.00 

24 

NOHINAL 

3 

IN 

0.00 

25 

(POSITIVE NHEN LARGER 

4 

IN 

0.00 

26 

THAN NOHINAL) 

5 

IN 

0.00 

27 

WSHATCH IN HHEEL ROLLING 

6 

IN 

0.00 

28 

AXLE 1 

IN 

0.000 

29 

RADIUS FROH SIDE TO SIDE 

2 

IN 

0.000 

30 

(POSITIVE NHEN HHEEL ON 

3 

IN 

0.000 

31 

LEFT RAIL IS LARGER THAN 

4 

IN 

0.000 

32 

THAT ON RIGHT RAIL) 

5 

IN 

0.000 

33 


6 

IN 

0.000 

34 

TRACK 

HALF KINEHATIC GAUGE 


IN 

29.75 

35 

RAILHEAD CROHN RADIUS 


IN 

10.00 

36 


EXHIBIT 3.3. 2.1 ~ VEHICLE PARAMETER OUTPUT 

(page 1 of 3 ) 
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ORiCINAl PASS is 

OF POOR QUALITY 


mss PROPERTIES 

LOCOMOTIVE BOOT 
HEIGHT 

ROLL HOHENT Of INERTIA 
PITCH HOHENT Of INERTIA 
YAH HOHENT OF INERTIA 


TRUCK FRAHES 
HEIGHT 

ROLL HOHENT OF INERTIA 
PITCH HOHENT OF INERTIA 
YAH HOHENT OF INERTIA 


HHEaSET/TRACnON HOTOR ASSEMBLY 
HEIGHT 


ROLL HOHENT OF INERTIA 
ABOUT CG 


OFFSET OF CG FROH 
AXLE CENTERLINE 

(POSITIVE HHEN HOTOR 
TRAILS AXLE) 

YAH HOHENT OF INERTIA 


LB 

ZS5700. 

37 

LB-IN-SECZ 

1510000. 

38 

LB-IN-SECZ 35300000. 

39 

LB-IN-SECZ 

35300000. 

40 


LB 

15440. 

41 

LB-IN-SECZ 

5Z655. 

4Z 

LB-IN-SECZ 

161400. 

43 

LB-IN-SECZ 

161400. 

44 


AXLE 1 

LB 

I31Z4. 

45 

Z 

LB 

131Z4. 

46 

3 

LB 

13124, 

47 

4 

LB 

131Z4. 

48 

5 

LB 

131Z4. 

49 

6 

LB 

13124. 

50 

AXLE 1 

LB-IN-SECZ 

ZOOOO. 

51 

2 

LB-IN-SECZ 

20000. 

5Z 

3 

LB-IN-SECZ 

20000. 

53 

4 

LB-IN-SECZ 

20000. 

54 

5 

LB-IN-SECZ 

20000. 

55 

6 

LB-IN-SECZ 

ZOOOO. 

56 

AXLE 1 

IN 

10.0 

57 

Z 

IN 

10.0 

SB 

3 

IN 

10.0 

59 

4 

IN 

-10.0 

60 

5 

IN 

-10.0 


6 

IN 

-10.0 

h 

AXLE 1 

U-IN-SECZ 

16780.0 

63 

n 

4. 

LB-IN-SECZ 

16780.0 

64 

3 

LB-IN-SECZ 

16780.0 

65 

4 

LB-IN-SECZ 

16780.0 

66 

5 

LB-IN-SECZ 

16780.0 

67 

6 

LB-IN-SECZ 

167B0.0 

68 


exhibit 3, 3. 2. I (CONTINUED) 
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SUSPENSION PfiOPERTIES 


OniGlNAi PAGE IS 

ccrnunADV OP QUAUtY 

SECONDARY 

LATERAL * 

STIFFNESS . *« 

RUBBER STIFFNESS 
SPRING TRAVEL TO STOPlEACH SIDE) 

VISCOUS COEFFICIENT (RUBBER) 

FRICTION BOLSTER COEFFICIENT 
LVB COEFFICIENT FOR BOLSTER FRICTION 
YAH DAMPING 

ROTATIONAL FRICTION FACTOR 
LVB COEFFICENT 
EXTERNAL DAMPERS 
VISCOUS COEFFICIENT 
LIMITING FORCE 


LB/IN 

15500. 

69 

IN 

1.250 

70 

LB/IN 

150000. 

71 

LB/( IN/SEC) 

136. 

72 


.40 

73 

LB/(1N/SEC) 

4000. 

74 

IN-LB/LB 

1.10 

75 

IN-SEC 

10000000. 

76 

LB/(IN/8EC) 

1020. 

77 

LI 

3600. 

78 


PRIMARY 

LATERAL 

STIFFNESS 

ROLLER BEARING FREE PLAY 
RUBBER ELEMENTS TRAVEL 
RUBBER PRELOAD 
STIFFNESS BEYOND STOP 
DAMPING 

FRICTION FOR ROLLER BEARING 
LVB COEFFICIENT FOR FRICTION ^ : 

VERTICAL 

STIFFNESS 

AVAILABLE TRAVEL(COMPRESSION) 
AVAILABLE TRAVEL(EXTENSIDN) 
SPRING~STIFFNESS(PER JOURNAL) 
STIFFNESS BEYOND STOP(EITHER HAY) 
PEDESTAL DAMPING 

FRICTION COEFFICIENTIPEDESTAL FACE) 
FRICTION COEFFICIENT(SIDE LUG) 

LVB COEFFICIENT FOR TOTAL FRICTION 


EXTERNAL DAMP.ERS (OPTIONAL) 
CODE TO INDICATE THE 
PRESENCE OR ABSENCE OF 
EXTERNAL DAMPERS 

(O^NO DAMPER) 
(IcDAMPER) 


RIGHT 


VISCOUS COEFFICIENT (COMPRESSION) 
VISCOUS COEFFICIENT(EXTENSION) 
LIMITING FORCE(COMPRESSION) 
LIMITING FORCE(EXTENSION) 


IS 


■M ^ 


LEFT AXLE I 
2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 


LB 

1500. 

81 

IN 

100000. 

82 


.10 

83 

LB/( IN/SEC) 

4000. 

84 

IN 

2.25 

85 

IN 

2.25 

86 

LB/IN 

5630. 

87 

LB/IN 

200000. 

88 

...» 

.’35 

89 


.35 

90 

LB/( IN/SEC) 

3000. 

91 


0. 

92 


1. 

93 

— 

0. 

94 


0. 

95 

— 

1. 

96 

— 

0. 

97 

— 

0. 

98 

— 

1. 

99 


0. 

100 

— 

0. 

101 


1. 

102 


0. 

103 

LB/( IN/SEC) 

510. 

104 

LB/dN/SEC) 

510. 

105 

LB 

IBOO. 

106 

LB 

1600. 

107 


HHEEL/RAIL 

aANGE 

aANGEHAY CLEARANCEfEACH SIDE) 
STIFFNESS AFTER FLANGE CONTACT 

TREAD 

CREEP COEFFICIENT ADJUSTMENT FACTOR 
COULOMB FRICTION COEFFICIENT 


IN 

LB/IN 


.550 lOfl 
60000. 109 

.70 no 
.30 111 


EXHIBIT 3, 3. 2.1 (CONTINUED) 


(page 3 of 3) 
3 -3 0 


2) Track Defect Output 


gjchiblt 3. 3. 2. 2 Is a sample output of the track 
defects specified for a particular run. In this case, 
three single defects, each 78’ long, have been defined. 

The first is a crosslevel defect which tilts to the left 
by Ij" in 59i"; the second is a lateral defect to the 
left of li", and the last is a vertical dip in the track 
of 1 3/8". 

3) Operational Parameter Output 

Exhibit 3. 3. 2. 3 is a sample output of the operational 
parameters . 
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NO. 

TYPE 

CLASS 

START 

AMPLITUDE HAVELENGH COUNT 




(FTI 

(IN) 

(FT) 

1 

X-LEVEL 

VERSINE 

100 . 

-1.250 

3S. 1 

2 

LATERAL 

VERSINE 

100. 

1.500 

3S. 1 

3 

VERTICAL 

VERSINE 

100. 

-1.375 

39. 1 


EXHIBIT 3. 3, 2. 2 - I’RACK DEFECT OUTPUT 

OPERATIONAL PARAHETER SPECIFICATION 


PARAMETER NAME 

VALUE 

UNITS 

LOCOMOTIVE SPEED 

40.00 

MPH 

NOMINAL TRACK CURVATURE 

3.00 

DEO 

NOMINAL TRACK SUPERELEVATION 

6.00 

IN 

TRACTIVE EFFORT ^ 

^9000. 

LB 

LATERAL HIND LOAD H 

a 0. 

LB 

LATERAL COUPLER FORCE/LEAD “ 

“ 0. 

LB 

LATERAL COUPLER FORCE/TRAIL 

0. 

LB 

RUN TIME 

5.0000 

SEC 

TINE STEP 

.0010 

SEC 

PRINT INTERVAL 

99.0000 

SEC 


EXHIBIT 3. 3. 2. 3 - OPERATIONAL PARAMETER OUTPUT 


FLAGS 0. 0. 0. 0. 0. 


TRUCK SIDE 
L/V RATIOS 

RATIOS TINES 
L R L R 



1 

.70 

.30 

2.12 

2.23 

.14 

.13 

2.24 

1.83 

.23 .21 

2.35 

2.37 


2 

.36 

.29 

2.40 

2.39 

.05 

.22 

2.50 

2.24 


3 

.26 

.16 

2.57 

2.59 

.07 

.07 

2.55 

2.99 





4 

.75 

.30 

2.B1 

3.01 

.14 

.12 

2.93 

2.63 

.22 .21 

3.04 

3.06 


5 

.38 

.29 

3.13 

3.07 

.05 

.20 

3.16 

2.90 



6 

.23 

.17 

3.29 

3.29 

.03 

.05 

3.27 

3.63 





EXHIBIT 3. 3. 2. 4 - SAFETY CRITERIA OUTPUT 


EUALUATION OF SAFETY CRITERIA 
PEAK VALUE ANALYSIS 


INDIVIDUAL NHEEL 
L/V RATIOS 

WTIOS TIMES 
R L R 


HHEELSET 
L/V RATIOS 


RATIOS 
L R 
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omGiNAt 

or POOR 

4) Summary of Results 

Exhibit 3 ^3. 2. 4 shows a sample summary output which 
Is divided Into two parts. 

fe. 

The first line, labelled "Flags", indicates abnovraal 
conditions which may have occurred during integration. 

The first flag is a count of Instances of wheel unloading, 
and the second is a count of instances of primary vertical 
unloading; the other three flags are not presently in use. 
Since the model assumes that the wheelsets follow the track 
and that net prlj^ary sus^nsion forces on each axle box 
are always downwards, remits from any run with non-zero 
counts are suspect. ^ 

i 

The bulk of the summary output reports peak values of 
the measures of performance, the various L/V ratios. For 
each of the L/V ratios, a maximum value and time is specified, 
For example, the peak L/V ratio that occurred on the right 
wheel of the second axle was .29 and occurred 2.39 seconds 
into the run. The peak wheelset L/V ratio on the fifth 
wheelset tending to shift the track to the left was .05 
at 3.16 seconds and to the right was .20 at 2.90 seconds. 


pfm 13 
quality 
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5) Output Outing mcegifa'tlion 

During integration there are two types of output which 
occur as described in section 3«1.4i these are generated 
by the routines SOUTPT and SPLOT. 

ft) SOUTPT 

Exhibit 3. 3, 2. 5*1 is a sample output from SOUTPT which 
is used for more detailed analysis, 

The first line indicates the time of the output and the 
current state of the five £iag^ discussed in above, 

fS ' t 

The second group provides instantaneous wheel/ rail 
forces as follows: 



FF 

hateral flange forces 

(lb) 


FCX 

Longitudinal Tread Creep forces 

(lb) 


FCY 

Lateral Tread Creep forces 

(lb) 


FWZ 

Vertical Wheel loads 

(lb) 


The 

twelve values given are for the six 

left wheels 

(axles 1 

through 6) and th- six right wheels 

(axles 1 through 

6) 

in that order. 



The 

next group provides, in the same order across 

the 

page: 




FPX 

Longitudinal primary fiorees 

(lb) 


Fpy 

Lateral primary forces 

(lb) 


FPZ 

Vertical primary forces 

(lb) 
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and these are follovrcd by the forces and moments at the 
secondary suspension for the lead and trail trucks: 


FSY 

Lateral secondary force 

(lb) 

FSZ 

Vertical secondary force 

<lb) 

MSX 

Secondary roll moment 

(in-lb) 

MSY 

Secondary pitch moment 

(in-lb) 

MSZ 

Secondary yaw moment 

(in-lb) 


Following these suspension forces and moments are 
listed the relative displacements (or rotations) and 
velocities at the interconnections: 


YW 

Wheelse t/rail lateral displacement 

(in) 

DYW 

Wheelset/rail lateral velocity 

(in/sec) 

THETW 

Wheelse t/rail yaw rotation 

(radians) 

DTHETW 

Wheelset/rail yaw velocity 

(radians /sec) 

YP 

Primary lateral displacement 

(in) 

DYP 

Primary lateral velocity 

(in/sec) 

ZP 

Primary vertical displacement 

(in) 

DZP 

Primary vertical velocity 

(In/sec) 

YS 

Secondary lateral displacement 

(in) 

DYS 

Secondary lateral velocity 

(in/sec) 


These are followed by the non-dimensional slip velocities 
of the wheelsets (WN) and the secondary yaw data: 


TRETS Secondary yaw rotation (radians) 

DTHETS Secondary yaw velocity ( radians / sec ) 
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The next three lines of data provide the ins tantaneoua 
data on the state variables (displacements and rotations, 
velocities, and accelerations)* 

The balance of the data relates to the L/V ratios: 

FL net lateral force outward on rail (lb) 

FV net vertical force on rail (lb) 

EC(12) individual wheel L/V ratios (ND)* 

EC(6) wheelset L/V ratios (ND) 

EC(4) truckslde L/V ratios (ND) 

* Non Dimensional 
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b) ai’LOT 


Xl\Q output to disc for off*-lino proousalng is not 
usually displayed, but must be roCerenced by index 
number for the off-line plotting. Exhibit 3>3. 2.5.1 
shows the index numbers. Exhibit 3. 3. 2. 5. 2 is a 
aatnplo output showing a plot of the body vertical 
displacement (2). As shown, the index for this variable 
is 181. 
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PlOniNO fUOCRAR 

RUN i: ! date: bo/ob/ii. 13.52.oz. 


ORlGINAt PAGE IS 
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IDHNT.'BDDY V.IN. SYfiBOL < • INDEX: • 181 

ti • 

^2.200 -1.760 -1.320 -.EfiO -.MO -.000 ,M0 .BBO 1.320 1.760 2.200 



EXHIBIT 3.3, 2.5. 2 - SAWLE PLOT OUTPUT; 

BODY VERTICAL DISPLACEMENT (IN.) 
VERSUS TIME (S) 
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3 • ^ COMPUTER RESOURCE REQUIREMENTS 


The program has been run on a CDG Cyber 174 computer 
under control of the NOS operating system, The definition 
of the required resources In terms of memory and computing 
tltne relates specifically to this machine, 

Although the time of execution depends upon the 
compiler In use, the rule of thumb Is that a typical FORTRAN 
program will run approximately uix times faster on the 174 
than on an IBM 360, This may help to provide a guideline to 
those not familiar with CDC hardware. 

On the CDG Cyber 174, the program requires approximately 
13 seconds of computer time and 60K words of memory for 
compilation and loading using the standard FORTRAN compiler. 

Execution of the program requires 18K words of memory 
and approximately 11 seconds of computer time per second 
of simulation. Additional time of about 0.75 seconds per 
second Is required for each defect and about 0,1 seconds 
per output (from SPLOT) , 

Define: 

t^ = total time to be simulated (seconds) 

Nj » total number o£ defects specified 

Np = number of SPLOT outputs per simulated second 

Then we can approximate the CPU time (t ^ ) by; 

cpu' ^ 

‘cpu = ^ («• ♦ <•”> 

On the CDC Cyber 174, our experience for typical runs 
indicates a cost of between one and two dollars per 
simulated second, 
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APPENDIX 4 


TECHNICAL DESCRIPTION OF THE MODEL 
4.1 INTRODUCTION 

The methodology contains an analytic model of; a 
locomotive which consists of a vehicle resting on two 
3-axle trucks. Figure 4.1 is a side-elevation view of 
the model (rear truck not shown), and Figure 4.2 is the 
corresponding end-view. These diagrams are designed to 
highlight the interface regions between the major functional 
parts of the locomotive suspension system, namely, the 
secondary suspension between the vehicle body and the truck 
frame, the primary suspension between the truck frame and 
the wheelsets, and the wheel-rail interface . 

Also shown in Figures 4.1 and 4.2, are some of the 
Important dimensions that have an Influence on the dynamic 
behavior of the locomotive. The definitions of these and 
other symbols that are being used in the modelling can be 
found in the technical description that follows. It is 
sufficient here to discuss a few of them, in order to point 
out some fundamental features in the model. 

- The axles of a truck are equally spaced; hence the 
half truck wheelbase A. 

The center of gravity of the truck frame is assumed 
to be located vertically above the center axle. 
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The center o£ gravity oi the vehicle Is assumed to be 
located midway between the center of gravity of the two 
trucks; hence the longitudinal truck center distance, 

2L. 


Whenever coupling is to be consldeted (Whether fore 
or aft or both) coupler forces will be acting On the 
corresponding coUpler-pln(s) ; the pins are assumed to be 
located at a distance D longitudinally from the center of 


gravity of the vehicle, and at a distance h^ above the rail- 
head . 


The tractive effort T of t%e locoirtotive Is assumed 
to be contributed by equal amounts at each the 6-axles. 

Braking force is treated as a negative tractive effort . 
Within a single test run, the tractive effort is assumed to 
be constant in magnitude. 


The locomotive is assumed to travel at a vollcity V 
which is constant duflng any single test run. 

The track has hOmlnal and constant values for curvature 
(D^) , absolute superelevation 0, and kinematic gauge 2b. 

Irregularities in the track geometry are supplied as 
input parameters. The irregularities covered are 
lateral alignment, vertical alignment c rOSS-leyel , and 
track gauge. 


As a result of disturbances transmitted via the track 
or other external forces, the locomotive (vehicle, trucks, 
wheelsets) is set in different modes of vibrations, arid forces 
are generated at the various Interfaces. These interfaces 
or suspension forces are treated in Section 4.3. Track 


geometry defects are described in Section 4.4, as well as 
other forcing functions such as coupler and wind loads. 

This Is followed by a description of the measures of 
safety in Section 4.5, and the equations of motion in 
Section 4.6. 

4.2 COORDINATE SYSTEM AND DEGREES OF FREEDOM 

Figure' 4.3 shows schematically the plan view of the 
locomotive at equilibrium on an ideal track with constant 
radius of curvature R . Because of the track curvature, 
the center of gravity of the vehicle and the trucks (A, B, and 
C, respectively), are all offset from the track. At this 
position of equilibrium, the center of gravity of each truck is 
displaced towards the center of the curvature of the track 
by an amount 



The reference systems and the degree uf freedom of the 
vehicle and of the trucks are defined as follows. 

1 ) Vehicle 

The origin of the Cartesian reference system coincides 
with the center of gravity of the vehicle, with the x-axls 
pointing along the longitudinal axis of the vehicle, and with 
* the x-y plane horizontal. See Figure 4.1 and Figure 4.3. The 

z-axis points vertically upwards, while the y-axls points away 
, from the center of track curvature. 

Since constant train velocity is assumed, there are 
five permissible degrees of freedom associated with the 
vehicle body : 
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y - vehicle lateral 
% - vehicle vertical 

a - vehicle roll 
0 - vehicle pitch 

Y - vehicle yaw 

The rotational movemanta follow the right-hand corkscrew 
sign convention about the k, y and z axes* respectively 

2) Trucks and Wheelsets 

The reference system of the leading truck frame and its 
wheelsets (identified as 1* 2, and 3* starting from the leading 
position) has an origin coinciding with the truck frame center 
of gravity. 

There are two degrees of freedom associated with the leading 
truck, frame: 


Pj^ - truck lateral, oriented along the radius of 
curvature OB and away from the center 0, 

6^ - truck yaw, about the vertical axis through the 

center of gravity, and observing the right-hand 
rule as shown* 


The wheelsets have one degree of freedom each, namely, 


lateral 
to the 


displacements y^^, y and y^, all 
radius of curvature, and measured 


defined 
from the 


pare lie 1 
center 


of gravity. 


The reference system of the trailing truck frame and its 
wheelsets, A, 5, 6, ate defined in a manner similar to that of 
the leading truck. 
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It ahould be noted tliat truck raoCiona in the vertical, 
pitch, and roll directiona are aaauroed to be constrained to 
those of the locomotive body at the secondary suspension. 

This assumption is considered valid when the secondary 
suspension stiffness is significantly higher than the primary 
suspension stiffness in those modes. Another assumption is 
that the wheelsets are constrained in yaw to the truck frame, 
This assumption is considered valid provided the longitudinal 
pedestal clearances are small. 

4 , 3 DESCRIPTION OF THE SUSPENSION ELEMENTS 

4.3,1 Do flnitions and N oinencl a t ure System 

The locomotive suspension elements can be classified 
with respect to three interfaces as follows: 


Interface 

Upper Part 

Lower Part 

secondary suspension 
primary suspension 
wheel-rail 

vehicle 
truck frame 
wheel 

truck frame 
wheelset 
rail-head 
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The noiiienclacure syscem uied to deClae the interface 
forces is as follows, 

Secondary suspension. 


^SX, 

k 

.pa 

secondary 

F 

®^k 

mm 

secondary 


- 

secondary 


Sm 

secondary 



secondary 

”-k 

PP' 

secondary 


Notes: F “ Force, M ^ moment^ s •“ secondary, 

X, y and z denote direction, 
k ”* 1, 2 for leading and trailing truck 
respectively . 

Primary suspension: 


F 

primary 

longitudinal force 

F 

=“ primary 

lateral 

force 


■" primary 

vertical 

force 

Notes: 

F ““Force, p 

primary 


X, y and z denote direction, 

i 1, 6 denotes Wheelset identification, 

j •* 1 , 2 for left and right sides respectively. 
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Whcol-rail Intorfoae: 



ilk 


vortical wheel-rail force 


^cx *** longitudinal creep force 

“lid 

F lateral creep force 

^^i,j 

F “ lateral flange force, 

^i.d 

Notes: P ** Force, w wheel, 

c creep, f •* flange, 

X and y denote direction, 
l,j same as for primary. 

To facilitate visualizing how these forces act in the 
locomotive model, an exploded free-body diagram is shown 
in Figure 4.4. 


4,3.2 Secondary Suspensions Characteristics 


1) Secondar y Lateral 

The secondary suspension lateral forces at the leading 
and the trailing trucks are each made up of four components: 

a) a component due to the nonlinear stiffness of the rubber 
elements and the limiting stops, 

b) a second component due to linearized viscous damping 
of the rubber material. 
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a third component due to the frictional properties of 
bolster on the truck traction pads^ 


d) a fourth component due to external hydraulic dampers 
when specified. 


The nonlinear stiffness of the secondary suspension is 
modelled as a piecewise linear function of the relative 
lateral displacement y^|^ as shown Ih Figure 4,5. It is 
characterized by a slope. 


K <iue to the rubber elements, 

sy 

and a much steeper slope when the spring t^ravel is 

exceeded and the stop la reached. 


The Coulomb friction properties of the traction pads is 
modelled as a function of relative lateral velocity in the 
secondary suspension With a linearized Viscous band (LVB) 
followed by saturation. (See Figure 4.6, where the LVB 
has a slope of ^ saturation of • T'he value of 

saturation depends on the limiting friction due to the 
longitudinal secondary force F 

0 X • 


2 ) 


Secondary Yaw M 


s z 


The secondary yaw moment M at either truck is due to 

s z 

the rotational ftlctiort encountered at the center-plate of 
the frusk. A LVB model similar to that of Figure 4.6 is used for 
this friction which la a function of the relative rotational velocity 
between the truck fralme and the vehicle body. The saturation 
level is Itself a linear function of the vertical secondary 
suspension force F > 
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3) Secondary tongitudinal F 


sx 


The secondary longitudinal forces acting at the 
center-plate of the leading truck ^nd of the 

trailing truck together balance the sum of all 

external forces acting on the locomotive. Assuming no 
dynamics for the locomotive in the longitudinal direction, 

^sxl ^sx2 *" ^ (total tractive effort . 

As mentioned earlier, a braking effort is treated like 
a negative tractive effort. 


4) 


Secondary Vertical F 


sz 


Assuming no Independent vertical dynamics of the 

individual truck frames, the vertical force F acting at 

s z 

the secondary suspension of the locomotive is equal to the 

sum of all the primary vertical forces F on a given truck. 

p z 


5) 


Secondary Roll M 


sx 


The secondary roll moment acting on the leading 

truck is the result of the sum of roll moments acting 
on the truck frame due to the primary vertical forces and 
primary lateral forces. 


6 ) 


Secondary Pitch M 


SI 


The secondary pitch moment at each truck is found by 
considering the balance of pitch moments acting on the truck 
frame due to the action of primary vertical forces, and the 
tractive effort. 
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4.3.3 Primal* y apena I on rac ter is tics 


1) Primary Lateral 


The lateral force F generated at each journal box 

py 

is due to a nonlinear stiffness characteristic in 
addition to friction forces developed in the roller bearingt 
The primary lateral stiffness is based on element tests on 
the resilient thrust unit. The test data are simplified and 
represented by the curve OAIBC as shown in Figure 4.7, 

The first portion Of this cUrve» OA, approximates a freeplay 
region of magnitude dj^. This is followed by a polynomial 
stiffness AB , covering a Working range of permissible bearing 
movement d^. The Value P. at A represents a preload on the 

I n 

rubber unit of typically 1500 lb. Beyond this range, a hard stop 
is reached, exhibiting a very steep load-displacement 
relationship BC . 


The frictional component of the primary lateral force 
is modelled again by a LVB as depicted in Figure 4.6. The 
saturation level, due to a limiting frictional force, is 
based on the resultant of tbe primary vertical and the primary 
longitudinal forces acting at tbe bearing in question. 

A coefficient of friction of O.i is normally used. 


2) Primary Vertical F 

-P ^ 


The vertical force <F acting at the end of each axle 

)p * 

is made up ot three coihpo'nen ts ; 
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a) tt stiffness component due to Spring stiffness and a 
limiting stop, 

b) a component due to pedestal face and pedestal side 
lug friction, 

c) a component due to an external damper, when specified. 


I’he stiffness is characterized by a piecewise linear curve 

as shown in Figure 4.8. It consists of a central segment 

of slope representing the spring stiffness. It is 

terminated at both ends by a steeper slope K representing 

s 

the end of vertical travel of the spring. The suspension 
is preloaded by an amount Fj^^ , which nominally is equal 
to a share of the total sprung load of the locomotive. 

Since there are 12 journals per locomotive, 


F^ = (W + 2W^) / 12 
Ko V t ' 


where = weight of vehicle 

= weight of truck frame and bolster, 
t 


The pedestal friction is modelled by a LVB, where 

the saturation level is made up of contributions of friction 

due to the primary longitudinal force F and due to the 

P ^ 

primary lateral force F . Thus 

py 


saturation level F. = u 

do '^x 


where p^andp^ are the frictional coefficients acting at the 


F 

+ p 

F 

px 

y 

py 


surfaces of contact in question. 
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Wherever there Is an external vertical damper In the 
primary suspension, the nonlinear characteristics in Figure 
4.9 are used. It is simply an asymmetrical LVB, with damping 
coefficients in one direction of travel, and in the 
other. The saturation levels may also assume different 
values ( f 3 » ^ * 

NOTE; PRIMARY VERTICAL FORCES FOR THE E-8 LOCOMOTIVE TRUCK 

In the case of the E-8 locomotive, different equations 
are used to represent the truck. The E-8 truck features 
complete "wheel load equalization" (equalizer bars and 
primary springs located at an intermediate location between 
the wheel sets) and a secondary suspension which is softer 
than the primary. For this truck configuration, therefore, 
the parameters used to define the primary stiffness are 
obtained from considerations of equivalent stiffnesses of 
springs connected in series as follows: 

1 = 1 4- 1 


where, k^ = equivalent vertical primary stiffness 

Used in the model 

kp = measured primary suspension stiffness 
kg = measured secondary suspension stiffness 

To account for the load equalization feature (i.e. wheel loads 
are not affected by quasi-static differences in crosslevel 
or surface defects within the truck wheelbase) ,^the vertical 
primary stiffness forces are assumed to be equal within each 
side of each truck (l,e. there can only be four different 
values of primary Vertical spring force per locomotive). 
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I’hese spring forces are obtained from average track cross- 
level and vertical track geometry under each truck. 

The velocity dependent forces (i.o. pedestal friction) 
are treated in the same manner as for a non-equalized truck. 

3) Primary Longitudinnl P 

. P.^_ 

The primary longitudinal forces are obtained by consider- 
ing force and moment balance on each wheelset in a horizontal 
plane passing through the wheelset geometric center, Longitu- 
dinal pedestal forces result from longitudinal creep forces, 
tractive effort^ and inertial reaction forces resulting from 
accelerations of the wheelse t / trac Cion motor assembly in Che 
lateral and yaw directions. 

4.3.^ Wheel-rail Interface 


1) Vertical Wheel Loads F 

W 2 


The vertical wheel-rail forces are obtained by consider- 
ing force and moment balance on each wheelset in a transverse 
vertical plane. Vertical wheel loads result from primary 
vertical forces, creep and flange forces, and inertial reaction 
forces resulting from track accelerations in the vertical and 
cross-level directions. 
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2) Wheel-Rail Creep forces F , F 
' ^ cx cv 


A non-linear creep force fornmljation is used, which 
is based on the wheel-rail contact geometry, the creepages 
at the contact points, the wheel-rail Coulomb friction 
coefficient, and the instantaneous normal tread force on 
each wheel. The formulation is: 

'• ■{■ ■ 3 ’}"" 

where: e = total creepage 

p = wheel-rail friction 
f = creep coefficient 
as shown in Figure 4,10* 

The creep coefficients are assumed to be equal in the lateral 
and longitudinal directions. This value is however dependent on 
the instantaneous wheel-rail vertical force as follows: 

f = G tr ab 

where G = shear modulus (11.5 it 10^ lb ) 

In^ 

a, b = semi-axis of the wheel-rail contact area 
(classic Hertz contact theory). 

Having obtained the total creep force, the longitudinal and 


lateral components are, 

respectively ; 

II 

• 

c. 

cx c 

X 


G 

a* 

II 

u. 

£ 

cy c 

y 


e 
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e being the vector sum of the longitudinal and lateral 

crecpage components e and e , respectively. The creepage 

X y 

components are computed from simple geometric relationships, 
taking account of the following; 

crack curvature and gauge, 
locomotive speed, 

wheelset angle of attack and angular velocity, 
wheelset deflection from track centerline and 
lateral velocity, 

- wheel tread conicity, rolling radius, and tape 
size mismatch as specified for each wheelset, 
wheelset rotational slip (resulting from 
tractive effort). 

3) Lateral Flange Contact Forces F^ 


When the wheelset is displaced laterally, the flange 
of one wheel touches the rail after the flangeway clearance 
is traversed. Further displacement will result in increas- 
ing contact force which can be modelled as a linear spring 
force. The flange clearance may deviate from the nominal 
value according to the local track gage variation at that 
point as described in the following section. 

4.4 FORCING FUNCTIONS 

4.4.1 Tra ck G eom etr y 

In addition to the nominal geometry, consisting of 
curvature superelevation and standard track gage, track 


geometry irregularities are considered^ Track geometry 
irregularities are defined as a summation of defects 
which are superimposed at the time of computation. Each 
defect is defined by the following: 

- Identification number (1,2, 3, etc.) 
starting point along track 
type (cross-level, lateral, vertical, or 
gage variation) 
class (verslne, or others) 
characteristic length (e*g> wavelength) 
characteristic amplitude 
number of repeated cycles. 

For each individually specified defect, the track 
input is computed using an expression of the form: 

(e.g. for verslne) 

‘I*. - Q J 1 cos 2r ( x - xs ) I I 

2 Y L (xe-xs )J J 

where x = position along track 

xs, xe ■“ start and end position of the defect 
Q — amplitude 

The Instantaneous position of each wheel set Is computed 
from the current simulated time, the locomotive speed, and 
the wheelset spaclngs. 

4 .4.2 Other F orcing Functions 

The other factors that are included in the model as 
additional Input functions arc; 
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tractive effort 


wind load F 

w 

coupler loads 

superelevation deficiency force F 

S 


1) Tractive Effort 


The tractive effort T acts on the locomotive in the 

e 

longitudinal direction at the three suspension levels. 
Traction provided by the traction motors is a positive 
tractive effort, while braking is regarded as a negative 
tractive force. 

2) Wind Load 

W 


A lateral wind force acting on the vehicle is considered. 
The wind load Is assumed to act at Che specified center 
of pressure for the locomotive. 

3) Lateral Coupler Forces F^ 


Lateral components of coupler forces can be specified 
at each coupler (i.e. leading and trailing ends). Presently, 
only a constant level of force can be specified which is 
applied throughout the simulation. 

4) Superelevation Deficiency Fg 


For a given superelevation, if the locomotive speed V 
is greater or less than the balance speed, superelevation 
deficiency prevails in the form of lateral loads acting through 
the center of gravity of the vehicle (Fg^) > the truck, frames 
(F ), and each of the wlieelset and traction motion subassemblies 

o ^ 


(F )• This load Is calculaxod as follows: 

8 w 

-g = «• <-r - »%> 

where m = mass 

V ~ train speed 

- nominal superelevation angle 

R =i radius of track curvature 
c 

g =s gravitational constant 

4.5 MEASURES OF SAFETY 

Measures of safety are computed simply from the 
Instantaneous values of net lateral and vertical wheel- 
rail forces. Individual lateral and vertical wheel forces 
are obtained as follows; 



where the sign of the lateral wheel force is defined 
as positive when directed away from track centerline, 

The measures o£ safety are computed as follows: 

1 ) Wheel L/V Ratio 
"Wheel L/V" = ^ 

This is petforraed for each wheel for a subtotal of 
12 measures of safety. 






# 


2 ) 


Wl ieelset L / V Ratio 
"Wheelset 1/V" « ^ 


left - ti right 
left « V right 


where left and right refer to the position of the wheel, 
This is performed for each wheelset. Since wheelset L/V 
ratio requires to be monitored for both positive values 
(towards outer rail) and negative values (towards inner 
rail) this gives a subtotal of 12 measures of safety. 


3 ) Tr uck Side L/V Ratios 

L lead + L center -f L trail 


"Truck side L/V" 


V lead + V center + V trail 


where lead, center, trail refer to the wheelset position 
within a truck. This is performed for each truck and 
each side for a subtotal of four measures of safety. 


In summary, a grand total of 28 measures of safety 
are monitored for peak value analysis. 

The instantaneous values of measures of safety are 
processed through a numerical averaging filter before analysis, 
for peak values. A simple "running average" technique is used 
over a time duration x^. The value of x^ is currently set 
to 50 msec, for all measures of safety, but can be altered 
as specified. 
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4 . 6 EQUATIONS OF MOTION 


/( , 6 , i Nome nela Cure 

The following nomenclature is used in the equations 
of motions. For the definition of the subscripting system 
of forces and moments, refer to sections 4.3.1 and 4.4.2. 


m locomotive body mass 

V 

m^ locomotive body mass including truck frame masses 

m truck frame mass 

ro,^ ^.5 mass of wheelset/ traction motor assembly 
Icx,Icy,lcz effective mass moment of inertia of 
locomotive body about the center of 
gravity, in roll, pitch, and yaw 
respectively 

» effective mass moment of inirtia of truck assembly 
in yaw 

g gravity constant 

hi - . height of secondary suspension above rail 
hji height of body center of gravity above secondary 
suspension 

h , vertical distance between coupler and body center 
d 

of gravity 

h »-> vertical distance between center of pressure 
w c 

(wind load) and body center of gravity 
height of coupler above rail 
2b « track kinematic gauge 
2i » truck wheel base 
2L truck center distance 
2D ™ distance between coupler pins 

os offset of wheelset/ traction motor assembly 
center of gravity from axle centerline. 


Sp 
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/i,6.2 Equations 


Tho 15 equations o£ motion are: 


1 ) Vehicl e La t e r a 1 


f, 


% 


I 


y^J:^<(^syl ■’’^’sy 2 > + Fao) » 


dl ‘d2' 


2 ) Veh i cle VettlGal Motion 

it 

2 ca ^ 




^szl ^ ^sz2? ” ® 


■]- 


3 ) Vehicle Roll Motion 
a 1 


( 


gv 


t «sx2> t »2 ‘‘'syl + >'ay2> 


+ "d <^dl + "d2) - V “'wind: 


A ) Vehicle Pitch Motion 


cy 


(“syl t ”,y2> t <^.2 




- T (h 
e p 


- h^) . 


5) Vehicle Yaw Motion 


^ [(M 

-cz { 


t M „) t L (F , - F 


■' 'T: ^ '"szl ’ “a^2’ ' “ ''syX " ‘sy2) 


+ ” <"dl - >’d2> 


] 
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wind 

] 


6) Leading Truck Lateral Motion 

3 2 

2 


- F 


m 


syl 


t js»l 


7 ) T railing Tjgu c k Lateral M ot ion 
r 6 2 

^2 ^ ’^pyij ' ’" ^'sy2 

™ t j“l 


8) Leading Truck Yaw Motion 

3 


’ 1 1 


Itf,! ^^cxl2 ~ ^cxll 


2 

^ ^ J^l ^^pylj 


M 


szl 


3 

23 

1«1 


2 

t : 

j»=«i 


^^cyij + ^fij 


F -, 4 ) 

py3j' 


%yij^ 


9) Trailing Truck Ya w Motion 

6 

0 lib ^ fF - F 

^ 1^4 '•^cxl2 cxll 


t2 


^ ^ j=.l ^'^pyAj ^py6j^ 


- M „ . c r s 


sz2 t 


^1 


j=Jl ^^cylj ^flj ^pyij^ 


]] 


b h 

10 to 15) 1 Wheelset/Tractlon Motor Subassembly Lateral 

Motion , 1 == 1 to 6 

2 


y e ® ^ 7 F 

^1 ®1 k i gwl 


m 


t jfl C^oylJ t - Fpyij 


] 


where k=>=l If l = x, 2, 3 
k== 2 If 1 =4, 5, 6. 
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• 6 . 3 Method of Solution 


From considerations of simplicity and convenience during 
program development and debugging, a simple Euler integra- 
tion technique was used. Determination of an adequate 
integration step size was performed empirically by sample 
problem tests over a wide range of step size values. An 
integration step size of .001 second was chosen as having 
a large margin of safety against perturbation of the 
resul t s i 
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FIGURE 4.1 -SIDE ELEVATION VIEW OF MODEL 
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SECONDARY SUSPENSION 








FIGURE 4.6 - LINEAR VISCOUS BAND (LVB) MODEL FOR 
COULOMB FRICTION 
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ORIGINAL PACE IS 



FIGURE 4.7 " PRIMARY LATERAL SUSPENSION CHARACTERISTICS 



FIGURE 4.8 


PRIMARY VERTICAL SUSPENSION CHARACTERISTICS 
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FIGURE 4,9 -EXTERNAL DAMPER CHARACTERISTIC IN THE 
VERTICAL PRIMARY SUSPENSION 
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FIGURE 4.10 - WHEEL-RAIL FRICTION CREEP CHARACTERISTIC 


-30 


